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PURPOSE 

The  basic  goal  of  the  total  study  program  is  to  establish  the  validity, 
capabilities  and  limitations  of  optical  power  spectrum  measurement  for  the 
determination  of  visual  display  information  capacity.  The  purpose  of  this 
initial  phase  is  to  establish  the  validity  of  the  basic  concept  and  ap- 
proaches to  be  used. 

APPROACH 

The  basic  measures,  transformations  and  analyses  were  examined  both  mathe- 
matically and  experimentally  with  respect  to  the  conditions  and  assumptions 
of  information  theoretic  measures  of  capacity.  The  basic  measurement  capa- 
bilities of  existing  equipment  were  assessed  in  terms  of  accuracy  and  repeat- 
ability. Potential  problems  in  practical  application  were  identified  and 
evaluated.  Possible  solutions  were  identified  and  evaluated  as  candidates 
for  further  development. 

Display  performance  parameters  were  reviewed  with  respect  to  effects  on  the 
information  capacity  measure  in  order  to  identify  measurement  requirements 
and  potential  areas  of  application  for  display  evaluation. 

RESULTS  AND  CONCLUSIONS 

The  general  approach  of  optical  power  spectrum  measurement  for  the  deter- 
mination of  display  information  capacity  is  valid  at  the  basic  level. 

The  basic  measurement  capabilities  of  the  existing  optical  power  spectrum 
measurement  equipment  and  calibration  procedures  are  adequate.  Root-mean- 
square  (RMS)  deviations  from  theoretical  values  are  about  10%.  RMS  measures 
of  repeatability  are  about  2%  of  the  mean.  These  values  apply  over  a useful 
dynamic  range  of  6 orders  of  magnitude. 
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There  are  a number  of  potential  problems  in  the  practical  application  of 
the  general  approach: 

o Nonlinearity  of  amplitude  units 
o Measuring  aperture  effects 
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o Phase  effects 
o CRT  film  recording 
o Noise  measurement 

Promising  solutions  for  these  problems  exist,  but  most  need  further  develop- 
ment and  testing. 

The  development  of  the  display  relationships  include  both  display  and  viewer 
parameters. 

o Spatial  frequency  response 
o Dynamic  range 
o Luminance 

o Visual  integration  time 
o Input  image  characteristics 
o Noise 

o Viewing  distance 


PREFACE 


This  report  was  prepared  by  the  Crew  Systems  organization  of  the  Guidance 
and  Control  Engineering  portion  of  the  Research  and  Engineering  Division 
of  The  Boeing  Aerospace  Company,  Seattle,  Washington.  The  work  was  done 
under  USAF  Contract  F33615-76-C-0030,  for  the  Visual  Display  Systems  Branch, 
Human  Engineering  Division,  of  the  Aerospace  Medical  Research  Laboratory. 

Mr.  Wayne  Martin  was  the  contract  monitor  for  AMRL. 

This  phase  of  the  study  was  started  in  May  1976  and  completed  in  July  1976. 
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SECTION  I 
INTRODUCTION 


The  Optical  Power  Spectrum  (OPS)  is  a physical  measure  that  can  be  used  to 
quantify  certain  characteristics  of  display  imagery.  It  is  an  attractive 
approach  because  of  its  speed  and  simplicity  and  because  it  operates  in  the 
frequency  domain  like  its  cousin,  the  modulation  transfer  function.  The 
latter  advantage  is  especially  important  when  working  with  displays  whose 
performance  is  frequency  dependent.  This  is  the  same  for  nearly  all  pictorial 
displays  such  as  photographs  or  television  receivers. 

The  OPS  measure  can  be  related  to  such  concepts  as  display  information  content 
and  capacity  as  developed  by  information  theorists.  Information  theory  has, 
for  many  years,  proved  to  be  a valuable  tool  for  the  analysis  and  evaluation 
of  comnuni cat ions  systems.  Application  to  pictorial  displays  has  been  develop- 
ing in  recent  years  and  promises  to  be  equally  useful  in  this  area.  The  pur- 
pose of  most  displays  is  to  provide  information  to  the  viewer.  Any  display 
is  limited,  spatially  and  temporally,  in  the  amount  of  information  that  it 
can  provide.  This  limit,  information  capacity,  is  an  important  parameter  in 
the  evaluation  of  display  quality.  Knowledge  of  the  impact  of  display  design 
parameters  on  information  capacity  and  of  the  impact  of  information  capacity 
on  viewer  performance  is  important  for  effective  development  and/or  selection 
of  displays.  The  measurement  of  information  capacity  is  a useful  tool  for  the 
development  of  this  knowledge. 

The  objective  of  the  total  study  effort  is  to  define  and  validate  OPS  measure- 
ment techniques  to  determine  display  information  capacity  and  to  apply  these 
techniques  to  the  evaluation  of  selected  display  parameter  effects  on  viewer 
performance.  This  report  covers  the  initial  phase  of  the  study  and  is  con- 
cerned with  the  validity  of  the  basic  OPS  approach  to  information  capacity 
measurement.  The  basic  approach  is  examined  both  theoretically  and  experi- 
mentally to  determine  the  capabilities  and  limitations  for  capacity  measure- 
ment and  to  outline  potential  areas  of  application.  Phase  II  will  develop 
the  specific  measurement  procedures  and  apply  these  to  the  preliminary 
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testing  of  the  relationship  between  the  capacity  measures  and  viewer  per- 
formance. Emphasis  will  be  on  CRT  displays.  Phase  III  will  apply  the 
resulting  procedures  for  the  evaluation  of  specific  display  parameter 
effects. 
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SECTION  II 

DISPLAY  INFORMATION  CAPACITY 


DEFINITION 

The  information  theory  approach  has  provided  a useful  tool  for  evaluating 
the  effectiveness  and  efficiency  of  information  transmission  systems. 
Application  of  the  theory  to  visual  displays  has  received  considerable 
attention  in  recent  years.  Reference  1 provides  a good  review  of  this 
work.  Appendix  A of  this  report  provides  a brief  background  of  informa- 
tion theory  and  the  development  of  the  basic  equation  for  display  informa- 
tion capacity. 

N log9L 

r _ ^ /i\ 


where  C = capacity  (bits  per  unit  display  area) 

N = number  of  independent  display  elements 
L = number  of  element  response  levels 
D = display  area. 

Equation  (1),  however,  is  limited  in  its  application.  Few  displays  of  in- 
terest have  independent  display  elements  (e.g.,  LED  and  liquid  crystal  dis- 
plays). The  response  of  a given  element  on  a CRT  or  photograph,  for  example, 
is  influenced  by  the  response  le.els  of  neighboring  elements.  The  number  of 
element  response  levels  is  usually  not  a fixed  value  but  often  varies  in  com- 
plicated ways.  A matrix  of  individually  illuminated  elements  can  meet  the 
assumptions  of  equation  (1),  i.e.,  independent  elements  and  fixed  number  of 
response  levels.  Even  here,  however,  the  assumptions  may  be  violated  with 
respect  to  the  viewer  if  the  viewing  distance  is  great  enough  or  if  the  illu- 
mination levels  are  at  viewing  threshold. 

A more  appropriate  definition  of  "N"  results  if  every  display  of  interest  is 
defined  as  "band-  imited."  That  is,  there  exists  a spatial  frequency,  K,  be- 
yond which  no  information  is  transmitted  to  the  viewer.  This  limit  may  be  im- 
posed by  the  display  (e.g.,  element  spacing)  or  by  the  capabilities  of  the 
observer's  visual  system  (i.e.,  visual  acuity).  Then,  the  minimum  number 
of  samples  needed  to  unambiguously  specify  any  given  display  response  is. 
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The  factor  of  4 is  necessary  because  of  the  requirement  for  2 samples  per 
cycle  in  each  direction.  Equation  (2)  applies  to  continuous  displays  such 
as  photographs  or  along  scan  lines  in  CRT  displays.  Equation  (2)  also  applies 
to  the  discrete  element  case,  where,  if 

d = spacing  between  elements 
then  K = 

The  factor  of  2 in  the  denominator  is  necessary  because  of  the  requirement 
for  2 samples  per  cycle.  Determination  of  the  number  of  response  levels, 

L,  is  the  most  difficult  and  controversial  aspect  of  informatirn  capacity 
measurement.  In  the  simplest  cases,  L is  determined  by  the  characteristics 
of  the  display.  This  is  true  tor  an  LED  matrix  display  where  the  elements 
are  either  on  or  off  and  the  illuminated  elements  are  easily  perceived  by  the 
viewer.  Then  L = 2.  In  more  complex  situations,  factors  such  as  noise  and 
system  frequency  response  characteristics  become  important.  If  the  display 
can  provide  any  output  level  over  a range  R . to  Rmax  then  L can  be  defined 
as 


L = |(R  - R . )/  AR 1 + 1 (3) 

max  min  I ' ' 

where  AR  = just  di scrimi nable  difference  in  signal 

response  level  R. 

Equation  (3)  is  in  the  form  of  the  familiar  signal-to-noise  ratio.  It  is 
strictly  valid,  however,  only  if  the  system  is  linear,  if  AR  is  independent 
of  the  value  of  R and  if  R„  , R . and  AR  are  constant  for  all  possible  dis- 
play  responses.  The  latter  requirement  is  generally  not  met  by  most  optical 
or  electronic  displays.  These  displays  are  usually  frequency  dependent  in 
their  response  levels.  Highly  detailed  (high  frequency)  patterns  have  a re- 
stricted range  of  allowable  response  levels  compared  to  gross  (low  frequency) 
patterns.  This  situation  requires  that  capacity  be  expressed  as  the  integral 


"f  the  frequency  related  L values,  so 

< 

x 


K 


C = 4 


J j log2  [L(k„,  Ky)] 


dkx  dky. 


(4) 


where 


kx  and  ky  are  spatial  frequencies  in  the  x and  y directions 
respectively,  and  K and  K , (the  upper  limits  of  integra- 

* y . 

tion)  are  the  band-limited  cutoff  frequencies  in  the  x and 
y directions. 


Determination  of  the  values  L (kx,  ky)  is  simplified  if  the  display  response 
function,  R (x,  y)  is  transformed  to  the  frequency  domain  by  using  the  Fourier 
transform. 


F(k 


V 


= //R(x,y  ) exp  [- 


T!  1 


(kxx 


kyy)]  d*  dy-  (5) 


Using  the  general  form  of  equation  (3)  we  can  define 


where 


L(kx,  ky)  = 


FHAX(t«’  V + , 


(6) 


FMAX^kx’  k ^ 1S  the  Four^er  transfonT1  of  the  maximum 
display  response.  FAR(kx>  k ) is  the  Fourier  transform 
of  the  just  discriminate  display  signal  response. 


Thus,  we  can  rewrite  equation  (4), 


loa 


■2 


FMAX(kx’  V 


dk  dk 
x y 


(7) 
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DETERMINATION  OF  ^AX ( kx , k ) 

Information  theory  demonstrates  that  F(k  , k ) is  maximized  when  the  input 

x y 

signal  for  the  display  has  equi-probable  element  response  levels  with  zero 
memory  (see  Appendix  A).  This  will  occur  when  the  input  signal  is  random 
over  a response  range  that  is  at  least  as  great  as  the  dynamic  range  capa- 
bilities of  the  display.  F^(kx>  ky),  then,  can  be  obtained  by  providing 
the  appropriate  random  input  to  the  display  and  determining  the  Fourier  trans- 
form of  the  resulting  display  response. 

It  may  also  be  meaningful  to  determine  capacity  for  selected  classes  of  input 
signals  such  as  alphanumeric  symbols  or  aerial  photography.  It  is  generally 
better  in  such  instances  to  use  a measure  defined  by  information  theory  as  "infor- 
mation content." 


H = CD 


where 


Then, 


H = information  content  (bits) 

C = capacity  (bits/unit  area) 

D = display  area 

h'4d <8> 


where 


F (k  , k ) is  the  Fourier  transform  of  the  displayed  signal 
s x y 

of  interest. 


DETERMINATION  OF  F^R(kx,  k ) 

The  most  common  approach  here  is  to  use  the  display  noise  level  as  the  basis 
for  measurement  (References  1,  3,  4).  The  response  distribution  for  a 

constant  level  input,  R^(x,  y)  is  Fourier  transformed  to  produce  F^(kx,ky) 
and  the  assumption  made  that 


F AR(kx’  V FN^kx’  V 
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and,  so 


iov 


FMAX*kx’  V 
W ky> 


* 1 


This  approach  assumes  that  performance  is  display  noise  limited.  It  is 
strictly  valid  only  for  systems  that  are  linear  in  amplitude  and  for  noise 
distributions  that  are  signal  independent  and  additive.  The  effects  of  de- 
viations from  these  conditions  are  not  well  known. 


There  are  other  approaches  available.  F D(k  , k ) could  be  determined  ex- 

AK  x y 

perimentally  for  a given  display.  A psychometric  evaluation  or,  perhaps, 
existing  visual  performance  data  could  be  used  to  establish  just  detectable 
display  response  levels  R^R(x,  y)  as  a basis  for  the  desired  values.  This 
psychometric  approach,  however,  is  time-consuming  and  expensive  and  the  re- 
sults limited  to  the  specific  displays  tested.  The  use  of  existing  data  is 
questionable  if  significant  display  noise  levels  are  involved. 

A compromise  approach  is  to  examine  F^(kx,  k ) and  its  effect  on  viewer  per- 
formance for  display  conditions  of  interest.  The  objective  of  such  an  approach 
would  be  to  determine  if  some  functional  form  of  the  transforms  would  be  a 
more  effective  descriptor  of  display  performance,  i.e.. 


f 

FMAX(kx’  ky). 

+ 1 

g 

TM(k  , k ) ‘ 
N x y' 

. 

j 

An  important  tool  for  this  approach  will  be  based  on  the  following  multiple 
regression  model, 


FINPUT^kx’  V = A + 8 [FS^x’  V]  + C [W  V]  • 

Cross  spectral  correlation  will  be  used  to  evaluate  the  functional  form  and 
appropriate  constants  for  this  model.  (See  Reference  5 for  a similar  approach. 
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The  resulting  correlational  statistics,  coefficient  of  determination  and 
standard  error  of  estimate  show  promise  as  indices  of  display  information 
loss  and  hence  as  indirect  measures  of  display  information  capacity.  This 
approach  will  be  reated  in  greater  detail  in  Section  V. 


DYNAMIC  DISPLAYS 


The  determination  of  F,(k  , k ) and  FM(k  , k ) for  dynamic  displays  such  as 
a CRT  presents  a special  problem.  The  noise  response  in  such  instances  is  prop- 
erly described  as  RN(x,  y,  t)  where  t is  a time  dimension.  Furthermore,  al- 
though we  are  concerned  here  only  with  fixed  input  signals,  the  display  re- 
sponse includes  the  effect  of  noise  and  will  vary  with  time.  Therefore,  we 
must  consider  R<.(x,  y,  t)  as  well.  The  effect  of  the  dynamic  noise  compon- 
ent will  be  modified  by  the  visual  integration  of  the  viewer.  If  T is  the 
visual  integration  time,  then 


F(kx,  ky,  T)  = 


j j P(x>  T)  exp  £-2  iri(kx  x+  kyy)J 


{kx  X+  kyy)|  dxdy  ^ 


X Y 


where 


R(x 


l 

, y,  T)  = J R(x,  y,  t)  dt 


t = 0 


and 


C = 4 


If 


log. 


FS,  n^v  V T2 
~F  AR^V  kv’  ^ 


dkx  dky 


(10) 


Note  that  in  the  above  expression,  the  numerator  F<-  ^ ( kx , k^ , T ) includes  the 
effects  of  noise  and  the  additive  factor  of  1 is  eliminated.  Thus, 


FS,N(kx,ky,T)  FS(kx’ky,T) 

^vVT> = F^Wr)  + 1 
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SECTION  III 

OPTICAL  POWER  SPECTRUM 


DEFINITION 

The  approaches  for  display  information  capacity  outlined  in  the  previous  sec- 
tion all  require  the  determination  of  the  Fourier  transform  of  a display  re- 
sponse pattern.  This  study  effort  is  concerned  with  the  application  of  a 
specific  technique  for  such  a determination.  The  optical  power  spectrum  is 
a physical  measurement  based  on  the  principles  of  coherent  optical  process- 
ing. Under  proper  conditions,  an  image  illuminated  with  coherent  light  pro- 
duces a diffraction  pattern  having  an  intensity  distribution  directly  related 
to  the  Fourier  transform  of  the  input  image  amplitude  transmittance  distri- 
bution. Specifically  under  the  conditions  shown  in  Figure  1, 


ry)  , 2C2  I FT^kx’  ky ’ I 


where,  I ( r , r ) is  the  diffraction  pattern  intensity  distribution, 
x y 

r , r are  linear  dimensions  in  the  transform  plane, 

x y 

L is  the  amplitude  of  the  coherent  illumination, 
x is  the  wavelength  of  the  coherent  illumination, 

F is  the  focal  ^ength  of  the  transform  lens, 

F.pj( kx , T)  | is  the  modulus  squared  of  the  Fourier  transform 
of  the  input  amplitude  transmittance  distribution  (after  integration 
over  time  T). 

A description  of  the  physical  process  and  the  development  of  equation  (11)  is 
presented  in  Appendix  B. 


The  optical  power  spectrum  P(k  , k , T)  is  defined  as, 

x y 


P(kx,  ky,  T)  = 


F(kx,  k , T) 
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1 


where,  P(kx,  kv,  T)  is  in  units  of  power  per  unit  frequency  souared 


D is  the  area  of  the  input. 


Combining  equations  (11)  and  (12) 


2 2 

p<>v  V T>  = ts:  ! (r*'VT) 


L D 


(13) 


It  is  important  to  note  that  the  optical  power  spectrum  (and  the  resulting 
capacity  measure)  is  based  on  the  amplitude  distribution  in  the  input  plane. 
This  fact  has  several  significant  consequences.  First,  the  amplitude 
distribution  of  the  display  image  is  probably  not  the  best  description  of  the 
display  response,  R (x,  y,  T).  The  visual  system  responds  to  intensity  (the 
square  of  amplitude)  not  amplitude.  Most  visual  performance  data  are  based 
on  intensity,  not  amplitude  units.  Furthermore,  we  can  expect  linearity  prob- 
lems, although  the  same  difficulties  will  exist  with  intensity  since  the  visual 
response  is  not  linear  with  intensity.  The  effects  of  this  problem  (amplitude 
vs.  intensity)  will  be  a major  concern  throughout  this  study  effort. 


A second  consequence  of  amplitude  dependence  is  that  amplitude  is  sensitive  to 
phase  effects.  Intensity  and  the  human  visual  system  are  not.  When  film  is 
used  as  the  input  medium,  phase  shifts  may  occur  because  of  variation  in  the 
film  thickness.  (References  6,  7,  8,  9,  10).  These  effects,  if  significant, 
must  be  eliminated  or  minimized. 


The  third  problem  is  that  the  power  spectrum  reflects  the  distribution  of 
the  total  input  plane.  This  includes  the  effect  of  any  aperture  used  to 
limit  the  image  area  being  measured.  The  Fourier  transform  of  the  limit- 
ing aperture  is  convolved  with  the  transform  of  the  input  image  and  a simple 
correction  is  not  feasible  (Reference  9). 


It  is  often  convenient  to  work  with  polar,  rather  than  linear  coordinates, 
in  the  diffraction  plane.  With  the  optical  axis  as  the  origin,  we  define 
the  angle  y and  the  distance  r as  shown  in  Figure  2.  Then 
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rx  = r cos  y 
ry  = r sin  Y 


and  (11)  becomes. 


Ur,  Y,  T)  = I F(k  , y , T)  I 

IT  1 r 1 


Equation  (12)  becomes, 

p(k  , Y,  T)  = 1 F(kr’  y’  T)  1 . 


So 


D 


2 2 

P(k  > Y.  T)  = I ( r , y , T). 

L D 


(14) 


(15) 

(16) 


The  value  of  1(0,  0)  is  of  special  concern.  This  value  is  related  to  the 


"DC"  or  average 

illuminance  value  taken  over  the  input  plane. 
1(0,  0)  = 1 F (0,  0)  T)  |2  . 

From  (11) 

F(0,  0,  T) 

= XfAT(x>  y.  T)  exP  ^"2  tt i ( 0 + 0)]dxdy 

= //at(x,  y,  T)  dxdy 
= 7^T  D 

where 

Ay  = the  average  amplitude  transmission 
D = input  image  area 

and 

L27Ld2 

1(0,  0)  = . 
kF^ 

(17) 

From  (13) 

2 2 

P(0,  0,  T)  = 1(0,  0) 

i/d 

= I2  D 

(18) 
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MEASUREMENT  REQUIREMENTS 


Equation  (16)  defines  the  basic  OPS  measurement  process.  The  intensity  of  the 

diffraction  pattern  as  a function  of  position  must  be  measured  and  then  mul- 
tiplied by  a constant  based  on  the  measurement  conditions. 

We  can  separate  the  required  measurement  system  into  several  functional  areas: 

o Illumination 
o Image  control 
o Optical  transformation 
o Diffraction  Pattern  Measurement 

Optical  considerations  and  the  analytical  conclusions  derived  above  define  a 

number  of  requirements  for  each  of  the  functional  areas. 

Illumination 

1.  The  input  image  must  be  illuminated  with  a coherent,  collimated  source. 

2.  The  level  of  illumination  must  be  uniform  within  the  system  limiting  aperture. 

3.  The  level  of  illumination  must  be  consistent  with  the  sensitivity  range  of 
input  image  transmission  values. 

4.  The  level  of  illumination  must  be  stable  over  reasonable  measurement  periods. 

Image  Control 

1.  The  input  image  must  exist  in  a form  capable  of  varying  the  illumination 
amplitude  or  phase.  Only  amplitude  variation  will  be  used  here.  Photo- 
graphic film  will  be  the  primary  medium.  Unwanted  phase  effects  must  be 
controlled. 

2.  The  image  must  be  positioned  so  that  the  system  limiting  aperture  covers 
the  image  area  of  interest. 

3.  The  image  position  must  be  stable  throughout  the  measurement  period. 

4.  The  image  should  be  positioned  perDendicular  to  the  optical  axis 

of  the  measurement  system.  This  is  not  a firm  requirement.  Deviations 
will  shift  the  diffraction  pattern  axes  in  a predictable  manner  but  opti- 
cal realignment  is  required. 
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Optical  Transformation 


1.  The  transform  lens  should  not  induce  significant  aberrations.  The  lens 
transfer  function  should  not  show  a significant  drop  in  response  over  the 
spatial  frequency  range  of  interest.  Transfer  effects  can  be  treated 
analytical ly; but  reduced  response  levels  create  signal-to-noise  problems 
in  the  diffraction  pattern  measurement. 

2.  The  transform  lens  focal  length  should  be  consistent  with  the  physical  size 
constraints  of  the  system  and  with  the  capabilities  of  the  diffraction  pat- 
tern measurement  components.  The  transform  lens  focal  length  determines 
the  "scale"  of  the  diffraction  pattern. 

3.  The  transform  lens  diameter  and  position  must  not  introduce  vignetting  over 
the  range  of  image  areas  and  diffraction  angles  of  interest. 

Diffraction  Pattern  Measurement 

1.  Diffraction  pattern  intensity  must  be  measured  at  the  transfer  lens  focal 
plane  and  relatable  to  position  in  that  plane. 

2.  The  integrating  area  of  the  individual  sensing  elements  must  be  sufficiently 
large  to  provide  adequate  sensitivity  and  small  enough  to  provide  reasonable 
spatial  frequency  resolution  in  the  resulting  spectra. 

3.  The  dynamic  range  capabilities  of  the  measurement  system  must  be  adequate 
for  the  expected  intensity  variation  in  the  pattern  (up  to  5 or  6 orders  of 
magnitude). 

4.  Accuracy  and  stability  of  the  intensity  measurements  must  be  consistent  with 
acceptable  error  levels  for  the  resulting  spectra. 

In  addition  to  the  specific  requirements  above,  good  optical  alignment  practices 

must  be  followed  to  insure  proper  geometry  and  distribution  of  the  diffraction 

pattern  and  to  minimize  the  effects  of  back  reflections  and  stray  light. 


EQUIPMENT 


All  optical  power  spectrum  measures  in  this  effort  will  be  made  with  a Record- 
ing Optical  Spectrum  Analyzer  (ROSA)  manufactured  by  Recognition  Systems,  Inc. 
This  equipment,  shown  in  Figure  3,  is  a general  purpose  instrument  designed  for 
ease  and  flexibility  of  optical  set-up.  The  optical  components,  all  with  mag- 
netic mounts,  are  positioned  on  a 72  x 44  inch  optical  table.  A fiberglass 
shroud  covers  the  table  and  components  to  protect  against  stray  light. 

For  this  work,  the  equipment  is  configured  generally  as  shown  in  Figure  4. 

The  illumination  source  is  a 7 mw  He-Ne  laser.  Illumination  is  controlled 
with  a variable  attenuator.  The  attenuator  consists  of  two  opposing  glass 
wedges  that  are  positioned  for  the  desired  transmission  without  imposing  a 
significant  deviation  in  the  laser  light  path.  The  laser  beam  is  expanded 
with  a microscope  objective.  The  magnification  of  the  objective  controls 
the  degree  of  expansion  which,  in  turn,  is  determined  by  the  desired  area 
of  illumination  at  the  input  plane.  A 25  pm  pinhole  is  placed  at  the  focal 
point  of  the  microscope  objective  to  "clean  up"  the  laser  beam.  The  laser 
and  microscope  objective  can  be  considered  as  a crude  spectrum  analyzer.  The 
pinhole  is  placed  at  the  Fraunhofer  diffraction  plane  of  the  objective  and 
thus  passes  only  the  uniform  energy  in  the  laser  beam.  It  acts  as  a spatial 
filter  that  blocks  the  high  frequency  variation  in  the  beam.  The  expanding 
beam  is  then  collimated  with  a lens  placed  at  a distance  from  the  pinhole 
equal  to  its  focal  length.  A 15-inch  focal  length  telescope  objective  is  gen- 
erally used  for  collimation.  During  alignment,  collimation  is  assured  by 
inserting  an  optical  flat  in  the  beam  "downstream"  from  the  collimating  lens 
and  adjusting  the  lens  focus  until  the  resulting  interference  fringes  are 
minimized. 

From  the  collimating  lens,  the  beam  passes  through  the  film  holder.  This 
holder  is  designed  to  accept  roll  film  up  to  9 inches  in  width  as  well  as 
slides  or  film  "chips."  It  is  equipped  with  a variable  aperture  to  control 
the  area  illuminated.  A series  of  discrete  circular  apertures  with  diameters 
from  1/8  inch  to  2 and  1/2  inches  is  available.  The  holder  includes  a pro- 
jection unit  for  viewing  the  illuminated  image.  This  unit  is  used  to  aid  in 
positioning  the  image  in  the  laser  beam. 
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Figure  3:  Recording  Optica I Spectrum  Analyzer 
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Figure  4:  Optical  Configuration 
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The  light  diffracted  by  the  input  image  is  collected  by  the  transform  lens 
and  focused  on  the  photodetector.  The  transform  lenses  used  here  are  simple 
spherical  telescope  objectives.  The  transform  lens  focal  length  determines 
the  distance  between  the  lens  and  the  detector  and  hence  the  "scale"  of  the 
diffraction  pattern. 

The  ROSA  uses  a 64  element  photodetector  array  for  measurement  of  the  dif- 
fraction pattern  intensity  values.  Each  element  is  a diffused  silicon  photo- 
detector on  a common  substrate.  The  array  as  diagrammed  in  Figure  5 consists 
of  32  "rings"  and  32  "wedges."  The  diffraction  pattern  is  symmetric  so  the 
array  is  split  essentially  in  half,  each  half  measuring  a different  "character- 
istic" of  the  pattern.  The  rings  measure  intensity  as  a function  of  frequency 
integrated  over  direction  and  the  wedges  measure  intensity  as  a function  of 
direction  integrated  over  frequency.  The  central  ring  (Ring  1)  is  a complete 
circle  and  measures  the  pattern  "DC"  level.  The  succeeding  rings  are  nearly 
half  circles  and  provide  measures  of  intensity  integrated  over  a frequency  band. 
Thus 

*1  "i 

I (Ring  i)  = JJ I(r,  y)  dr  dy  (19) 

o m . 

l 

where  m. , n^  = minimum  and  maximum  radii  of  ring  i respectively 
f.  = angle  covered  by  Ring  i 

r,  y are  polar  coordinates  in  the  diffraction  plane. 

The  lowest  and  highest  frequencies,  a.  and  b^ , covered  by  Ring  i are  deter- 
mined from  equation 

ai  = m./xF 
b.  = n./xF 
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where 


X = laser  wavelength 
F = transform  lens  focal  length 


Values  for  iik  , n^,  y.  and  total  ring  area  are  listed  in  Table  1.  A 0.8  mil 
wide  conductor  strip  is  incorporated  in  each  ring  to  improve  detector  per- 
formance. The  effect  of  this  strip  on  total  ring  effective  area  is  included 
in  the  tabled  values. 

Table  2 shows  the  corresponding  spatial  frequency  values  for  the  40  inch  trans- 
form lens  used  in  this  study. 
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Table  1:  Detector  Ring  Dimensions 


Ring  No. 


Minimum 
Radius  (M. ) 
(mils)  1 


Maximum 
Radius  (N- ) 
(mils)  1 


108.8 


139.3 

157.6 

178.2 

201.2 
226.9 

255.7 

287.7 
323.2 
362.6 

406.0 

453.8 

506.4 

564.0 


107.8 


138.3 

156.6 

177.2 

200.2 
225.9 

254.7 

286.7 
322.2 
361.6 

405.0 

452.8 

505.4 

563.0 

626.0 


Ring  Angle 
(¥*) 

(degrees ) 

360.0 

180.0 
180.0 
180.0 
180.0 

166.3 
168.5 

170.1 

171 .4 

172.5 

173.4 

174.2 

174.9 

175.5 

176.0 

176.5 

176.9 

177.3 

177.6 

177.9 

178.1 

178.3 
178.5 

178.7 

178.8 
179.0 


179.2 

179.3 

179.3 

179.4 


Total  Effective 
Area  (A-j ) 
(sq.  mils) 


122.7 

148.6 

190.0 

245.0 

318.7 
417.4 

548.7 

722.0 

948.7 

1279.6 
1688.8 

2237.0 

2963.6 

3884.7 

5114.3 

6697.1 

8715.4 

11209.0 
14385.5 

18473.4 

23448.5 

29630.1 

37361.6 

46614.1 
57988.0 
71889.4 

88448.7 
108409.8 
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Table  2:  Detector  Spatial  Frequency  Values 


40  inch  Focal  Length 

Ring  No. 

Mi nimum 
Frequency 
(cy/mm) 

Maximum 

Frequency 

(cy/mm) 

Center 

Frequency 

(cy/mm) 

Bandwidth 

(cy/mm) 

1 

0.00 

0.14 

0.07 

.14 

2 

0.23 

0.34 

0.28 

.11 

3 

0.38 

0.49 

0.43 

.11 

4 

0.53 

0.65 

0.59 

.12 

5 

0.69 

0.80 

0.75 

.11 

6 

0.85 

0.97 

0.91 

.12 

7 

1.01 

1.15 

1 .08 

.14 

8 

1.18 

1.33 

1.26 

.15 

9 

1.37 

1.51 

1.44 

.14 

10 

1.57 

1.75 

1.66 

.18 

11 

1.79 

2.00 

1 .89 

.21 

12 

2.03 

2.27 

2 15 

.24 

13 

2.31 

2.57 

2.44 

.26 

14 

2.61 

2.92 

2.76 

.31 

15 

2.96 

3.31 

3.13 

.35 

16 

3.35 

3.75 

3.55 

.40 

17 

3.79 

4.25 

4.02 

.46 

18 

4.29 

4.81 

4.55 

.52 

19 

4.85 

5.45 

5.15 

.60 

20 

5.49 

6.17 

5.83 

.68 

21 

6.22 

6.99 

6.60 

.77 

22 

7.03 

7.90 

7.46 

.87 

23 

7.93 

8.91 

8.42 

.98 

24 

8.95 

10.05 

9.49 

1.10 

25 

10.08 

11.31 

10.69 

1.23 

26 

11.34 

12.70 

12.02 

1.36 

27 

12.75 

14.26 

13.50 

1.51 

28 

14.30 

15.97 

15.13 

1.67 

29 

16.01 

17.85 

16.93 

1 .84 

30 

17.89 

19.93 

18.91 

1.96 

31 

19.97 

22.20 

21 .08 

2.23 

32 

22.24 

24.68 

23.46 

2.44 
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The  "wedges"  are  identical  except  for  the  two  middle  ones  (Wedges  16  and 
17  numbering  clockwise).  All  wedges  measure  over  an  included  angle  of 
5.625  degrees.  All  have  an  outer  radius  of  626  mils.  Wedges  16  and  17 
have  an  inner  radius  of  26  mils,  the  others  of  21.4  mils.  The  total  effec- 
tive area  of  wedges  16  and  17  is  1748.8  square  mils  and  that  of  the  other  wedges 
it  is  1795.9  square  mils. 

Intensity  values  for  the  wedges  are  given  by. 


I (Wedge  i ) 


y)  dr  dy. 


ai  mi 


(20) 


where,  a-,  8.,  are  minimum  and  maximum  angles  of  wedge  i respectively, 
itk , r.  are  minimum  and  maximum  radii  of  ring  i respectively. 

The  outputs  from  the  individual  elements  are  fed  into  a pre-amp/mul tiplexer , 
then  each  signal  is  fed  into  an  auto-ranging  amplifier  where  it  is  converted 
to  the  proper  range  for  the  binary-coded  decimal  converter.  The  signals,  con- 
verted into  3 digit  binary-coded  decimals  are  then  recorded  on  punched  cards 
by  a Univac  Model  1710  card  punch  and  interpreter.  The  output  values  repre- 
sent microvolts  with  a maximum  value  of  1 x 10^  y v . Each  card  contains  a 6 
digit  designator  and  the  output  values  for  8 elements.  Eight  cards 
are  produced  for  each  measurement.  Time  for  the  measurement  sequence,  after 
image  positioning,  is  paced  by  the  key  punch  and  takes  about  20  seconds.  The 
data  taking  procedure  is  controlled  with  the  entry  and  control  panel.  This 
unit  enables  the  manual  entry  of  the  6 digit  designator  code  for  recording  on 
the  punched  cards. 

Data  processing  for  converting  the  raw  voltages  to  the  appropriate  power  values 
is  accomplished  with  a Hewlitt  Packard  Model  9820  calculator.  The  calculator 
is  equipped  with  a card  reader  for  data  input,  an  expanded  tape  memory,  and  a 
plotter  and  typewriter  for  data  output. 
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DATA  PROCESSING 


The  specific  processing  techniques  to  be  used  depend  on  the  nature  of  the 

input  image.  If  the  input  variation  is  one-dimensional  such  as  a single 

slit  or  a line  grid,  most  of  the  diffracted  energy  is  confined  to  a single 

spike  and  the  power  spectrum  is  properly  expressed  as  one-dimensional  (power 

per  cycle/nirn).  It  cne  image  information  is  truly  two-dimensional,  then  pro- 

2 

cessing  must  yield  the  two-dimensional  power  spectrum  (power  per  (cycle/mm)  . 
One -Di men sional  Variation 

The  individual  ring  voltage  values  are  converted  to  relative  intensity  per 
mil  after  correction  for  amplifier  bias  errors.  Bias  error  values  are  deter- 
mined with  a dark  current  (covered  detector)  reading  and  subtracted  from  the 
element  voltages.  Relative  intensity  is  determined  by  dividing  by  the  ele- 
ment width  as  determined  from  Table  1.  Except  for  ring  1,  the  elements  must 
he  corrected  for  the  .8  mil  conductor  strip  along  each  ring.  So, 


V-  - B. 

MA'  = , 


(21) 


where  L(A)  = relative  intensity  per  mil  for  ring  i, 

Vi  = raw  voltage  output  for  ring  i 

Bi  = dark  current  value  (bias  correction)  for  ring  i, 
n^,m*  = maximum  and  minimum  radii,  respectively,  for  ring  i. 


Then,  except  for  ring  1, 

2a  F ( V • - B ) 

P(kr,  T)  • «FI.(A)  m—  a- rfcr- 


(22) 


where,  F - the  focal  length  of  the  transform  lens  (mils), 

A = the  wavelength  of  the  laser 
* 

(6.328  x 1(J  mm/cycle  for  the  He-Ne  laser), 
the  center  frc 
(see  Taole  2) , 


k^  = the  center  frequency  for  ring  i 
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The  value  2 is  included  since  the  rings  2 through  32  cover  only  half  of  the 
total  pattern.  Ring  1 is  a full  circle  and  does  not  contain  a conductive 
strip.  The  center  frequency  for  Ring  1 is  considered  to  be  zero.  Thus, 


*F(V.  - B.) 

P(0,  T)  = A F I-j  (A)  = — . (23) 

' i 


Equations  (21),  (22),  and  (23)  provide  relative  values  in  that  they  do  not 
include  the  effects  of  laser  amplitude  or  image  area  as  required  by  equation 
(16).  This  is  generally  solved  by  normalizing  to  P(0,  T). 


2 (V.  - B.) 

P(k  , T)  ItTTTX"  .8 

p'(kr’ T)  = piorrr = "(v^  - b^ 

_____ 

1 


(24) 


The  prime  designation  is  used,  hereafter,  to  designate  the  normalized  power 
values.  Normalization  by  the  ring  1 value,  although  a common  practice,  is  not 
always  desirable.  First,  the  normalization  eliminates  the  influence  of  the 
average  image  transmission  level.  This  effect  may  well  be  of  interest  in  some 
applications.  Second.  the  normalization  should,  theoretical ly,  be  made  to  the 
true  "DC"  level  or  the  intensity  at  zero  frequency.  Since  ring  1 is  of  finite 
area,  it  includes  the  effects  of  very  low  frequencies  as  well,  and  thus  is  only 
an  approximation  of  the  true  "DC"  level.  Selection  of  the  best  normalization 
technique  will  remain  a major  concern  throughout  this  study  effort. 


The  wedges  are  not  of  direct  value  in  the  one-dimensional  analysis  except  to  verify 
the  appropriateness  of  one-dimensional  processing.  The  wedge  at  the  spike  orien- 
tation should  be  significantly  greater  than  the  general  level. 

Two-Dimensional  Variation 

Processing  is  extended  to  the  two-dimensional  input  by  using  the  element  area 
and  squaring  the  >F  correction  factor. 
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P(kr,  T) 


X2F2(Vi  - Bi) 


(25) 


where  A.  is  the  element  area  as  listed  in  Table  1.  These  values  include  the 
compensation  for  the  conductive  strip  and  for  the  variation  in  angular  cover- 
age of  the  rings  so  the  factors  of  2 and  0.8  used  in  equation  (22)  are  not 
required  and  equation  (25)  applies  to  all  rings  including  ring  1. 


Normalization  to  ring  1 provides, 

(Vi  - Bi) 

, P(kr.  T)  A 

P^(kr’  T)  “ 7 (0,  T)  = V1  - Bn  ' 

A 


(26) 


hocessing  for  the  wedge  data  is  complicated  by  the  reduced  areas  in  wedges 
16  and  17.  Area  correction  is  not  appropriate  here  since  this  correction 
assumes  an  equal  energy  distribution  across  the  element.  This  is  usually  an 
acceptable  assumption  for  the  rings  in  the  two-dimensional  analysis  but  does  not 
hold  for  the  wedges.  Standard  processing  techniques  for  the  wedge  elements 
have  not  been  established. 

Visual  Angle  Transformations 

When  the  power  spectrum  is  to  be  related  to  visual  performance,  units  relating 
to  visual  angle  are  preferred.  This  is  accomplished  by  the  use  of  the  follow- 
ing relationship. 


k 


e 


arctan  (1/D)  ’ 


where, 


k,  = spatial  frequency  in  cycles/degree 

b 

D = viewing  distance  in  millimeters 
k^  = spatial  freauency  in  cycles/millimeter 


and, 
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(28) 


P'(k  , T)  - p'(kr>  T) 


Larctan  TVD)  ] 


o 

where  P'(k  , T)  is  in  units  of  relative  power  per  (cycl e/degree ) . 


Intermediate  Image  Scale 

The  optical  power  spectrum  measurement  requires  that  the  input  image  exist  as 
a transparency.  Measurement  of  CRT  display  imagery,  for  example,  will  require 
the  use  of  a photograph  of  the  screen.  The  scale  of  this  intermediate  image 
must  be  used  to  adjust  the  dimensions  of  the  measured  spectrum  to  refer  to  the 
original  display  dimensions. 


SECTION  IV 
EXPERIMENTAL  TESTS 


INTRODUCTION 

This  section  describes  a series  of  tests  designed  to  provide  experimental 
demonstrations  of  the  general  validity  of  the  basic  approaches  described 
in  the  previous  sections.  The  tests  are  not  intended  as  an  evaluation  of  any 
specific  technique  for  the  measurement  of  image  information  capacity.  Results 
can,  at  best,  indicate  the  potential  for  practical  application  of  the  approach 
described  in  this  report.  Evaluation  of  specific  practical  techniques  is  the 
subject  of  Phase  II  in  this  study  effort. 

BASIC  POWER  SPECTRUM  MEASURES 

The  validity  of  the  basic  diffraction  integral  (equation  11)  can  be  examined 
by  considering  simple  apertures  as  input  images.  This  is  a useful  starting 
point  since  a simple  clear  aperture  in  opaque  material  represents  the  total 
input  distribution.  There  are  no  additional  limiting  aperture  effects  and 
no  phase  effects.  Transmittance  levels  are  0 and  1 and  so  provide  the  special 

situation  where  amplitude  and  intensity  are  equal.  Comparison  of  measured  with 
theoretical  clear  aperture  values  provides  an  indication  of  the  magnitude  of 
the  errors  that  may  exist  in  the  practical  application  of  the  theory.  The 
principal  sources  of  error  here  are  in 

1)  The  quality  and  alignment  of  the  optics, 

2)  The  accuracy  and  sensitivity  of  the  photodetector  and  electronics, 
and 

3)  The  accuracy  and  effectiveness  of  the  data  processing  and  calibra- 
tion techniques. 

One  Dimensional  Variati on 

We  can  express  the  definition  of  the  power  spectrum  from  equation  (12)  in 
one  dimension  as. 
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P(k  , T)  =IF"V  T>l 


where 


F(kx,  T) 


J R(x,  T)  exp  [ -2tti' kxx  ] dx 

size  of  the  image  in  the  direction  of 
interest  (x). 


The  simplest  case  of  a one-dimensional  measurement  is  a single  long  slit. 
Consider  a clear  slit  of  width,  w,  placed  symmetrically  about  the  origin 
and  uniformly  illuminated  so  that  the  response  function 


r(x,  T)  = 1 for  I x | s «/2 

I 0 elsewhere 


Since  the  aperture  is  time  invariant,  the  T designation  will  be  dropped. 
Over  the  area  of  integration, 

R(x)  = 1 

/W2  , r n+ai/2 

exp  (-2irikxx)kx  = - [exP  (-2rri kxx)  J 

h)‘  ^ sin  ( 7rkxo)) 


/ sin(-nk  w ) 

P(k  ) = I . — 

x I ukx(o 


^sin(nkxio) 


P(k  ) l irk  fci 
P (kx)  = P jOj  = ^sTn~TY2~ 
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where  — is  defined  as  equal  to  1. 

Figure  6 is  a plot  of  equation  (29)  for  w = 1.5  mm.  Since  the  power  spectrum 
is  proportional  to  the  diffraction  pattern  intensity,  the  plot  shows  the  slit 
diffraction  pattern  as  a line  of  closely  spaced  spots  decreasing  in  intensity 
from  the  center.  With  the  exception  of  the  center,  the  spots  are  spaced  0.67 
cycles/mm  apart.  Since  distance  in  the  diffraction  plane,  r,  is, 

r = F:\kr, 

the  spots  will  be  about  0.017  inches  apart  in  the  pattern  produced  by  a 40 
inch  transform  lens. 

The  finite  width  of  the  detector  ring  elements  will  integrate  the  energy  over 
the  corresponding  distances.  The  detector  thus  performs  a "smoothing"  fur  Jon 
in  the  measurement  process.  Since  the  ring  widths  increase  with  distance,  the 
smoothing  is  greater  at  higher  frequencies.  Smoothing  has  the  danger  of  hiding 
real  variations  in  the  pattern, but  is  useful  in  averaging  out  random  dis- 
crepancies to  produce  a more  stable  and  representative  spectrum.  Equation  (29) 
has  been  integrated  over  the  appropriate  frequency  ranges  for  the  individual 
rings,  assuming  a 40  inch  transform  lens  focal  length,  and  the  results  shown 
in  Figure  7.  The  effects  of  the  smoothing  process  are  clearly  evident. 

Five  repeated  optical  power  spectra  measurements  were  made  using  a 1.5  mm  wide 
slit.  The  ROSA  was  configured  with  a 40  inch  focal  length  transform  lens.  Mean 
values  were  used  with  equation  (24)  to  provide  the  data  points  shown  in  Figure  8. 
The  solid  curve  represents  the  theoretical  values  constructed  by  joining  the 
midpoints  of  the  bars  in  Figure  7.  The  root  mean  square  error  is  about  11%  of 
the  mean  value  over  the  range  of  more  than  4 orders  of  magnitude.  The  theo- 
retical values  are  strictly  correct  only  for  a slit  that  is  infinitely  long. 

The  slit  used  here  was  slightly  over  11  mm  in  length.  The  length  does  con- 
tribute to  the  diffraction  pattern.  The  magnitude  of  this  contribution  was 
calculated  and  found  to  be  3%  or  less  and  has  been  neglected  here.  The  repeat- 
ability of  the  readings  was  assessed  by  calculating  the  standard  deviation  for 
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Theoretical  Spectrum  ( 1.5mm  Slit) 


each  of  the  individual  elements.  Values  were  generally  of  the  order  of  1% 
to  2%  of  the  mean  reading. 

Wedge  readings  confirm  the  adequacy  of  the  one-dimensional  measurement.  The 
wedge  corresponding  to  the  spike  produced  by  the  slit  length  was  3.2%  of  the 
wedge  corresponding  to  the  major  direction  of  measurement.  The  repeatabi 1 ity 
of  the  wedges  was  consistent  with  the  levels  observed  for  the  ring  elements. 

These  results  demonstrate  the  validity  of  the  basic  diffraction  integral  as 
expressed  in  equation  (11).  They  also  indicate  that  the  measurement  equip- 
ment is  operating  satisfactorily  and  that  the  data  processing  procedures  are 
adequate  for  this  application. 

Two  Dimensional  Variation 

Consider  a clear  circular  aperture  with  diameter,  d.  Then, 


l 1.  |X|,  (V I £ d/2 
Rl*-  v>  ’ jo,  elsewhere 

As  before,  the  T designation  will  be  dropped  since  the  aperture  is  fixed. 

The  Fourier  transform  of  this  function  is  aerived  in  most  optics  textbooks, 
(see  References  2 and  11),  and  will  not  be  repeated  here.  The  result  is  the 
well  known  Airy  pattern  distribution.  In  polar  coordinates  this  result  can 
be  stated  as 


F (kr. 


y)  = F (0) 


2J.|  (ndsina/x ) 
Tidsina/A 


where 


a = the  angle  of  diffraction 
d = the  diameter  of  the  aperture 
( ) = first  order  Bessel  function 
F(0)  = zero  frequency  or  "DC"  amplitude. 


(30) 
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From  (B.l)  in  Appendix  B, 


sin  a = krA 
and  (30)  can  be  rewritten  as. 


2J  Mk  ) 

F kr*  >>  - F<°> 

r 


From  (15), 


P(kr,  y) 


F(kr,  y) 


F(0) 


2J1  (ndkr) 
iidk 


where  D is  the  area  of  the  aperture. 


Therefore,  since. 


P(0)  = 


F ( 0 ) 


then, 


P'(kr,  y)  = 


P(kr,  Y 

p(oT 


) 


| 2J1  (7Tdkr) 

\ '73k  I 

r 


(31) 


Equation  (31)  was  integrated  over  the  areas  of  the  individual  elements  to  derive 
the  theoretically  expected  values  for  three  circular  apertures  1/4,  1/2,  and  3/4 
inches  in  diameter.  Five  repeated  OPS  measurements  were  made  of  each  aperture 
with  a 40  inch  focal  length  transform  lens.  The  mean  readings  were  processed  as 
shown  in  equation  (26)  and  the  results  shown  in  Figures  9,  10,  and  11  along  with 
the  theoretical  curves.  The  most  obvious  point  in  the  figures  is  the  large  de- 
viations that  occur  at  power  levels  below  1 x 10_f\  It  appears  that  the  dynamic 
range  capabilities  of  the  ROSA  will  not  permit  accurate  measurement  below  this 
limit.  Performance,  however,  is  satisfactory  over  6 orders  of  magnitude.  The 
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Figure  9:  X Inch 


Figure  10:  & Inch  Diameter  Circular  Aperture 
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root  mean  square  error  levels  for  values  above  1 x 106  are  about  10%  for  the 
1/4  inch  aperture  and  about  20%  for  1/2  and  3/4  inch  apertures.  The  very  steep 
gradients  in  the  diffraction  pattern  for  the  larger  apertures  present  unusually 
difficult  measurement  problems.  The  1/2  inch  diffraction  pattern  intensity, 
for  example,  drops  6 orders  of  magnitude  in  a distance  of  less  than  0.2  inches. 
Fortunately  gradients  of  this  magnitude  rarely  occur  in  practical  applications 
and  the  10%  error  level  should  be  representative  of  system  performance.  The 
standard  deviations  for  the  repeated  readings  were  consistent  with  those  ob- 
served with  the  slit. 

For  the  basic  power  spectrum  measurements  then,  in  both  one  and  two  dimensions, 
the  equipment  and  data  processing  techniques  can  be  expected  to  provide  root 
mean  square  error  levels  of  + 10%  and  a repeatability  error  cf  + 2%  over  a dynam 
range  of  6 orders  of  magnitude. 

FILM  INPUT  MEASUREMENTS 

The  previous  section  dealt  with  the  performance  of  power  spectrum  measurements 
of  clear  apertures.  Most  measurements  of  interest,  however,  will  utilize  photo- 
graphic film  as  the  input  medium.  Figure  12  presents  an  example  of  an  OPS  mea- 
surement of  a photograph  of  a line  grid.  The  grid,  ts  measured  on  the  photo- 
graph, had  a frequency  of  1.16  cycles/mm  with  a line  (clear  space)  width  cf 
0.37  mm.  The  one-dimensional  theoretical  spectrum  can  be  expressed  as,  (after 
Reference  2), 


2 2 
n/i  ^ ( sin(nookx)  ^ ! sin(NuSkx)  ^ 

P(kx}  ” \ ™kx  / \ si nT^Skxr  / 

where  w = width  of  the  clear  spaces 

S = width  of  one  cycle 
N = number  of  bars. 

Comparison  of  the  measured  and  theoretical  spectra  in  Figure  12  indicates  good 
agreement  in  shape,  at  least  for  the  first  five  harmormcs.  The  discrepancies. 
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however,  exceed  the  accuracy  levels  found  for  the  clear  aperture  measurements. 
The  film  imagery  introduces  several  additional  sources  of  error: 


1.  Nonuniformities  in  the  thickness  of  the  film  Dase  and  emulsion 
introduce  unwanted  phase  effects  that  distort  the  spectrum. 

2.  The  aperture  used  to  limit  the  measurement  area  in  the  film  has  an 
effect  on  the  total  power  spectrum. 

3.  Film  granularity  introduces  a noise  effect. 

4.  Geometric  distortion  in  the  film  introduces  dimensional  measurement 
errors. 

5.  The  transfer  characteristics  of  the  camera-film  combination  modify 
the  contrast,  generally  as  a f notion  of  frequency. 

These  errors  will  be  treated  in  greater  detail  in  Phase  II.  However,  a pre- 
liminary test  to  evaluate  the  significance  of  the  phase  effects  was  performed 
in  this  series. 

As  discussed  in  Section  III,  the  optical  power  spectrum  reflects  the  input  am- 
plitude distribution  and  as  such  is  sensitive  to  phase  effects  introduced  by 
thickness  irregularities.  There  are  two  sources  for  this  problem;  the  variations 
in  the  film  base  and  the  variations  in  the  emulsion  caused  by  the  developed 
image  (image  relief  effects).  The  latter  problem  is  the  most  serious  since  it 
is  "signal  dependent"  and  cannot  easily  be  measured  independently. 

Liquid  gating  is  generally  recommended  to  eliminate  this  problem  (References  6, 

7,  8,  9,  10).  This  technique  involves  immersion  of  the  film  sample  in  a liquid 
of  similar  refractive  index  contained  in  a glass  cell  with  optically  flat  sides. 
Although  there  are  some  difficulties  because  the  film  base  and  emulsion  generally 
have  different  refractive  indices,  the  technique  is  successful.  However,  it  is 
a cumbersome  and  time  consuming  procedure.  The  matching  of  refractive  in- 
dices must  be  done  carefully  and  separately  for  each  film  type  used.  The  errors 
introduced  by  the  liquid  and  the  cell  must  be  determined  and  removed.  Because 
of  these  problems  it  is  important  to  determine  the  magnitude  of  the  phase  effects 
and  consider  otier  approaches  to  compensate. 
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An  image  containing  a number  of  bars  of  various  sizes  and  orientations  was  mea- 
sured and  then  bleached  in  photographic  reducer  to  remove  all  of  the  developed 
silver  grains.  The  resulting  image,  although  clear,  retained  the  original 
thickness  variations.  The  image  relief  was  clearly  visible  under  reflected 
light  and  in  transmitted  coherent  light.  Figure  13  compares  the  power  spectra 
of  the  original  and  bleached  images.  Both  measurements  were  made  with  the  lar- 
gest possible  anerture  size,  i.e.,  3/4  inch  and  a 40  inch  focal  length  transform  lens. 
Figure  14  compares  the  bleached  image  spectrum  with  that  from  a sample  of  clear 
film.  Although  it  is  obvious  that  the  image  relief  does  contribute  to  the  power 
spectrum  the  difference  is  small  when  compared  to  original  image  power  spectrum. 

The  difference  between  the  bleached  image  spectrum  and  the  clear  image  spectrum 
is  about  3 % of  the  total  image  power  levels  except  at  very  low  frequencies 
(<  0.5  cy/mm)  where  considerable  spreading  of  the  DC  spot  occurs  with  the  clear 
image.  This  comparison  is  not  intended  to  be  a definitive  experiment,  but  it  does 
suggest  that  the  use  of  a clear  film  sample  for  phase  effects  correction  may  be 
sufficient  for  practical  levels  of  accuracy.  Such  corrections  are  to  be  explored 
and  evaluated  in  Phase  II  of  this  study. 

INFORMATION  CAPACITY 

The  power  spectrum  provides  the  measurements  required  for  the  information  capa- 
city measure  as  defined  in  Section  II,  equation  (10), 


C = 4 


Jf 


log. 


'V'v'y1' 


dkxdky 


From  (11)  in  Section  III 


P(kx,ky,T)  = 


F(kx,k  ,T) 


and. 
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Figure  13:  Phase 
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C = 4 


//log. 


ps.N(kx,kv,T) 

PAR(kx’VT^ 


dkxdky 


= 2 


PS,N^kx’ky ,T^ 
PAR  kx  ky 


dkxdk 


y 


where  the  upper  limit  of  integration  is  the  system  cut-off  frequency,  K. 
In  polar  coordinates. 


C = 2 


PS.N(kr,y-T) 

PAR(kr’Y,T^ 


dk  dy. 
r ' 


(32) 


In  order  to  evaluate  the  application  of  equation  (32)  a set  of  imagery  was  pre- 
pared using  a computer  controlled  laser  beam  recorder.  The  recorder  is  a drum 
facsimile  type  controlled  by  a Varian  6201  computer  which  is  responsible  for 
timing  and  logical  operations.  Digital  image  data  are  provided  via  magnetic  tape. 
The  computer  performs  a gamma  correction  which  compensates  for  nonlinearities 
in  the  printer  modulator  and  film.  The  data  are  then  converted  to  analog  form  to 
drive  the  laser  modulator. 


The  imagery  consists  of  nine  frames  with  random  dot  patterns.  The  dots  are  0.1 
mm  wide  with  0.1  mm  spacing.  The  frames  represent  three  different  numbers  of 
grey  levels  (2,  4 and  16)  each  at  three  levels  of  dynamic  range.  The  grey  levels 
are  printed. at  equal  density  intervals  and  are  based  on  equal  probabilities  of 
occurrence.  The  steps  are  thus  linear  in  density  but  have  a log  relationship 
with  intensity  and  amplitude.  All  nine  conditions  were  printed  to  the  same  den- 
sity level.  Table  3 lists  the  magnitude  of  the  density  steps  for  each  condition. 


No.  of  Levels 


Table  3: 


Grey  Step  Sire 

(Density  Units) 

Dynamic  Range 

High 

Medi urn 

Low 

0.90 

0.30 

0.15 

0.30 

0.10 

0.05 

0.06 

0.02 

0.01 

2 

4 

16 
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Table  4 defines  the  dynamic  range  conditions. 


Table  4:  Levels  of  Dynamic  Range 


Density 


Range 

Maximum 

Minimum 

Mean 

Range 

High 

1.25 

0.35 

0.80 

0.90 

Medium 

0.95 

0.65 

0.80 

0.30 

Low 

0.88 

0.72 

0.80 

0.16 

A tenth  frame  was  produced  with  a uniform  signal  at  the  same  mean  density  level. 

The  film  was  Type  2405  XX,  developed  in  D19  for  7 minutes  at  70°  F.  The  printer 
calibration  curves  show  a linear  response  over  the  maximum  density  range  used 
(0.35D  to  1.25D). 

Optical  spectrum  measurements  were  made  on  all  10  frames  using  a 40  inch  trans- 
form lens  focal  length  and  a 1 inch  diameter  circular  aperture.  The  resulting 
power  spectra  normalized  to  ring  1 are  shown  in  Figure  15.  Except  for  the  med- 
ium contrast-2  level  frame,  the  spectra  do  not  vary  with  the  number  of  grey 
levels  but  are  sensitive  to  dynamic  range.  The  spikes  occur  at  the  frequencies 
corresponding  to  the  spacing  of  the  dots.  The  first  spike  occurs  at  about  10 
cycles/mm  corresponding  to  the  inverse  of  the  dot  spacing  at  orientations  of  0° 
and  90°.  The  second  spike  is  at  14.1  cycles/mm  and  represents  the  second  harmonic 
of  the  diagonal  sampling  freauency. 

The  power  spectra  should  not  vary  with  the  number  of  grey  levels.  Power  levels 
are  a function  of  the  input  amplitude  variance  (mean  square  amplitude)  and  since 
the  dot  pattern  grey  levels  have  equal  probability  of  occurrence  the  resulting 
distributions  (rectangular)  all  have  equal  theoretical  variance.  The  discrepancy 
for  the  medium  contrast-2  level  frame  is  most  likely  attributable  to  physi- 
cal alignment  inaccuracies  in  the  laser  beam  recording  process. 

The  information  capacity  of  the  input  signals  can  be  calculated  from  equation  (1), 
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Figure  15:  Random  Level  Dot  Spectra 


where,  for 


N - lumber  of  dots  = = 100 

L - number  of  response  levels  {2,  4,  or  16) 
2 

D-d  red  - 1 mm  . 


Tab1?  b lists  these  results. 

Table  5:  Information  Capacity  of  Input  Siena Is 


Numbei  of  grey  levels 


Capacity  (bits/mm*-) 

100 

200 

400 


The  informal  .'or.  cc-rcity  of  the  display  (laser-film  combination)  can  be  de- 
termined as  di scir  :;y-  In  Section  it  by  using  the  uniform  density  power  soectruns 
as  a measure  .v  Then  for  a noise  limited  display,  we  ran  define  the 

)ust-disr»*i::iin«i rspc.nse  level  in  terms  of  the  noise  and  equation  (32)  be- 
comes. 


- i / • 

“ p*c  .1  ] 

^ * : Ova 

\j  d 

‘■w-r  | 

J 

(33) 


The  prime  desi-jual  ions  have  yeer,  added  to  indicate  the  use  of  norma  1 ved  power 
spectra.  Mcma;  spectra  are  not  required  if  all  measurements  are  made  with 
the  same  laser  em,:!  i-;»de  and  r.he  >a:ne  iiraq?  moon  transmission  values.  Although 
these  confi  or  ■ controlled . corua  "Sized  soectra  were  used  here  tc  correct 

for  possible  discrepancies.  Th;  T designation  (integration  time)  has  been  dropped 
since  this  analysis  evolves  a static  display. 


it  is  important  to  note  that  equation  (33)  defines  the  capacity  of  the  display. 
It  does  not  reflect  the  amount  of  information  actually  displayed.  As  with  the 
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power  spectrum,  this  measure  should  not  vary  with  the  number  of  grey  levels 
in  the  signal  but  should  be  a function  of  the  dynamic  range.  Capacity  values 
were  calculated  for  each  of  the  nine  frames  using  equation  (33).  The  results  are 


listed  in  Table  6.  The  value  K = 
the  upper  limit  of  integration. 


1 


2 x dot  spacing 


= 5 cycles/mm  was  used  as 


p 

Table  6:  Display  Capacity  Values  (bits/mm  ) 

Number  of  Input  Grey  Levels 

Dynamic  Range  2 4 16  Mean  S. D. 


0.90D 

282.1 

288.7 

289.1 

286.6 

3.9 

0.30D 

109.2 

138.6 

128.8 

125.5 

15.0 

0.16D 

54.2 

56.6 

53.2 

54.7 

1 .8 

The  agreement  across  input  levels  is  very  good  except  for  the  discrepant  frame. 
The  mean  capacity  values  vary  as  expected  with  dynamic  range.  If  the  dynamic 
range  is  expressed  as  modulation,  the  relationship  between  information  capacity 
and  dynamic  range  is  nearly  linear  over  the  values  measured. 


m = ^max  ~*min 

T + T . 
max  min 

where  T , T . are  the  maximum  and  minimum  intensity  transmission  values, 
max  min 

T = antilog  (-D). 

This  is  shown  in  Figure  16.  The  amount  of  information  for  each  of  the  three  in- 
put levels  is  shown  as  dashed  lines.  The  plotted  values  indicate  that  the  dis- 
play capacity  is  exceeded  by  the  16  level  input  at  all  levels  of  dynamic  range. 
Capacity  is  exceeded  by  the  4 level  input  at  the  medium  and  low  dynamic  range 
levels  and  by  the  2 level  input  by  the  low  dynamic  range  level. 

The  preceding  analysis  demonstrates  that  the  general  behavior  and  reliability  of  the 
OPS  information  capacity  measure  support  the  validity  of  the  basic  approach. 
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Consider  the  4 grey  level-high  dynamic  range  conditon.  From  Figure  16,  the 

p 

input  information  level  (200  bits/mm  ) is  within  the  display  capacity.  From 
Table  3,  the  grey  step  size  is  0.30D.  Also  from  Table  3,  the  2-level  medium 
range  condition  represents  a single  grey  step  for  the  former  condition. 

Thus, 


C = 2 


If  '°92 


ww> 

P2M(kr,Y.T) 


dk  dy , 
r 


should  equal  the  input  information  capacity,  where 


P4H(kr'f,T)  = power  spectrum  measured  on  the  4-level  high  dynamic 
range  condition. 

i°2M^V’y,'0  = power  spectrum  measured  on  the  2 -level  medium  dynamic 
range  condition. 


Evaluation  of  the  above  equation  provides  a value  for  the  output  image  informa- 

o 

tior  capacity  of  180  bits/mm  . This  value  compares  favorably  with  the  theore- 

p 

tical  value  of  200  bits/mm  . 
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SECTION  V 


IMAGE  NOISE 


INTRODUCTION 

The  use  of  the  optical  power  spectrum  for  determining  display  information 
capacity  is  highly  dependent  on  the  proper  measurement  of  image  noise.  When 
the  display  system  is  noise  limited,  the  display  noise  power  spectrum  is  a 
major  component  of  the  information  capacity  metric.  A major  requirement  is 
that  the  noise  measurement  be  truly  representative  of  the  form  and  magnitude 
of  the  noise  actually  present  in  the  image  of  interest.  The  ease  with  which 
this  requirement  can  be  met  is  largely  determined  by  the  characteristics  of 
the  noise  itself. 

Power  spectrum  measurement  of  photographic  noise  has  been  extensively  studied 
(see,  for  example,  ref.  6 and  10).  Photographic  noise,  principally  granularity, 
is  one  of  the  least  difficult  types  to  evaluate.  It  is  static,  isotropic 
(equal  in  all  directions),  behaves  in  a reasonably  predictable  fashion  and  has 
been  examined  in  great  detail.  This  is  in  contrast  to  CRT  display  ncise,  for 
example,  that  is  dynamic,  differs  significantly  with  direction,  and  is  not 
well  understood.  There  has  been  a great  deal  of  work  with  noise  in  the  CRT 
video  signal;  but  much  less  is  known  about  the  resulting  image  noise  on  the 
CRT  screen. 


The  following  discussion  considers  those  noise  characteristics  that  are  important 
for  the  proper  measurement  and  use  of  the  noise  power  spectrum  in  it. formation 
capacity  measures. 


ADDITIVITY 


The  information  capacity  equation  that  uses  the  noise  power  spectrum  as  a 
direct  signal -to-noise  ratio, 


c ■ ‘ // 


dk  dk 

x y 


n y 

is  based  on  the  assumption,  among  others,  that  the  noise  tehaves  in  an 


(34) 
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additive  fashion.  This  means  that  the  noise  is  added  to  the  signal  level. 
Thus , 


RS,  N {X’  y)  = RS  (X’  y)  + RN  (X’  y) 
and  from  the  properties  of  the  Fourier  transform 


Fs,»  (kx-  V ' FS  "V  V + F„  <kx'  ky> 


and  since 


P (kx.  ky)  = 


F ( k , k ) 
' x y' 


ops,n  (kx 


'V  |FS.N(kx'  ky>  I ■ lF5<kx'ky>  + [rN<Vk,>  I 


" |FS<kx’k»>l  +2Fs<kxV  FN.(kx'V  + lF»'kxV 

* 

• DPS<kx-k,)tDFn(kx-k,)t2Fs<kx'ky»FN(kx’ky> 

(35) 

If  the  noise  is  also  independent  of  the  signal,  the  cross-product  term 
(representing  covariance)  becomes  zero  and  the  power  spectra  are  also  additive. 

If  the  display  noise  is  multiplicative  rather  than  additive;  i.e., 


R (x,  y)  = R (x,  y)  R (x,  y) 
b >N  5 N 


then  it  is  the  case  that 


FS,N  (kx’  ky)  = lFs  (kx>  VI  lFN  ky) 


and 


PsVkx’  ky^ 


Fs  (kx>  V 


-12 

F (k  , k ) | 

N X’  y>  | 


(36) 


where  * denotes  the  operation  of  convolution. 
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Equation  (36)  indicates  a much  more  complex  effect  of  noise  on  the  power 
spectrum  and  the  use  of  equation  (34)  for  information  capacity  measurement 
is  questionable.  A theoretical  solution  for  multiplicative  noise  does  not 
exist.  If  the  errors  introduced  by  the  inappropriate  use  of  equation  (34) 
are  excessive,  then  modifications  in  the  form  of  the  equation  will  be 
required.  This  concern  is  not  limited  to  the  optical  power  spectrum  of  in- 
formation capacity;  it  applies  to  all  applications  of  the  signal-to-noise 
ratio  as  a description  of  system  performance. 


SIGNAL  DEPENDENCE 

One  aspect  of  the  correlation  between  signal  and  noise  has  been  mentioned  above. 
If  a correlation  exists,  the  cross  product  term  in  equation  (35)  does  not 
vanish  and  becomes  a source  of  error.  The  magnitude  of  the  error  depends  on 
the  size  of  the  correlation. 


Another,  although  related,  difficulty  occurs  with  the  selection  of  a representa- 
tive noise  sample  for  measurement  of  the  noise  power  spectrum.  If  the  noise 
is  truly  independent  of  the  signal,  then  the  noise  will  be  identical  at  all 
signal  levels  and  any  level  of  uniform  signal  input  can  be  used  to  obtain  an 
appropriate  noise  spectrum.  If,  however,  the  noise  level  varies  with  the 
signal  level,  the  decision  is  more  difficult.  The  common  solution  here  is  to 
use  theroise  measurement  at  the  mage  mean  signal  level.  This  may  not  be  the 
best  choice,  particularly  if  the  relationship  betv/een  signal  level  and  noise 
level  is  not  linear.  It  is  important,  therefore,  to  evaluate  the  magnitude 
and  the  form  of  the  relationship  and  select  the  most  appropriate  conditions 
for  the  noise  measurement. 

DIRECTIONALITY 


If  the  magnitude  or  other  characteristics  of  the  noise  vary  with  orientation 
in  the  display,  it  becomes  important  to  consider  this  effect  in  any  evaluation. 
There  are  two  approaches  to  this  problem.  One  is  to  use  only  input  signals 
that,  are  directionally  uniform.  This  will  avoid  the  interaction  between  scene 
content  and  noise.  The  other  approach  is  to  use  an  input  signal  with  highly 
directional  content  such  as  a grid  pattern  and  obtain  measurements  over  a range 
of  orientations.  The  latter  approach  will  provide  a quantitative  evaluation  of 
the  directionality  effect. 
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DYNAMIC  NOISE 


If  a display  output  varies  with  time,  the  effect  of  noise  that  is  random  with 
time  will  depend  on  the  magnitude  and  rate  cf  variation  and  on  the  integration 
time  of  the  viewer  or  other  sensor  (such  as  photograpnic  film).  If  the 
magnitude  of  the  noise  is  described  at  any  time  t by  its  variance,  then 


2 

a 


N,  t 


£(RX  - R)2 

n - 1 ' 


where 

Rx  = the  display  response  level  at  point  x 

R - the  average  response  level 

n = the  number  of  display  points  sampled 


If  the  noise  is  directional,  sampling  should  be  only  in  the  direction  of  interest. 


If  the  noise  is  random,  then  the  response  levels  at  each  point  in  trie  display 
will  be  averaged  over  some  integration  time  T,  In  general, 

2 2 

a a 

N,  T = N,  t 
n' 


where 


n'  = the  number  of  time  samples  in  the  integration  period. 


For  example,  if  f is  the  frame  time  for  a CRT  display,  then 


Parseval’s  theorem  states  that 

J |f(x)l  dx  = j |f(u)|  ciu 

For  this  application  the  above  is  interpretec  to  mean  that  the  area  under  the 
square  of  the  Fourier  transform  (i.e.,  power  spectrum)  is  equal  to  the  variance 
of  the  input.  Hence  a reduction  in  noise  variance  by  a tactor  n’  because  of 
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integration  will  produce  a reduction  by  the  same  factor  in  the  area  under 
the  power  spectrum.  It  may,  therefore,  be  possible  to  derive  the  spectrum 
of  the  integrated  input  analytically. 

The  empirical  approach  is  to  perform  the  integration  photographically  (by 
control  of  exposure  time)  and  perform  the  power  spectrum  measurement  on  the 
integrated  signal.  Both  approaches  are  to  be  considered  in  Phase  II. 

EVALUATION 

The  magnitude  of  the  effect  of  deviations  from  the  assumptions  of  additive 
and  independent  noise  is  difficult  to  evaluate  in  the  general  case. 

To  obtain  some  estimate  of  the  potential  error  magnitude,  a "worst  case" 
condition  was  tested. 

The  effect  of  overprinting  a "noise"  pattern  on  a photograph  is  to  produce 
multiplicative  noise  when  measured  as  transmission.  A test  of  such  noise 
effects  was  performed  using  photographic  transparencies  of  two  high  contrast 
line  grids  (1.16  cy/mm  and  1.76  cy/mm),  and  two  "noise"  overlays.  The  overlays 
were  produced  by  the  laser  recorder  with  a spot  size  and  spacing  of  0.2  mm. 

The  spots  varied  in  transmission  according  to  a random  Gaussian  distribution. 
The  mean  transmissions  of  the  two  overlays  were  equal.  They  differed  in  root- 
mean-square  deviation  (high  noise  and  low  noise).  The  dot  pattern  of  the 
overlays  was  aligned  with  the  bars  on  the  grid  to  produce  a high  spatial 
correlation  between  the  bar  signal  and  the  noise.  Optical  power  spectrum 
measurements  were  made  of  each  of  the  grids,  each  of  the  overlays,  and  of 
each  grid  under  both  overlay  conditions.  The  measurements  were  made  with  a 
40  inch  focal  length  lens  and  a 24  mm  x 35  mm  rectangular  aperture.  If  the 
noise  behaved  like  independent,  additive  noise,  then  in  polar  coordinates, 

PS,  N 'kr’  7 ^ “ PS  (kr*  ’ ^ + PN  ^kr’  7 ^ • 

For  each  of  the  grids  the  two-dimensional  power  spectra  (integrated  across  y ) 
were  calculated  and  used  to  evaluate  this  relationship.  Figure  17  illustrates 
this  evaluation  for  the  1.16  cy/mm  grid  and  high  noise  condition. 

The  spectra  were  evaluated  out  to  the  cutoff  frequency,  K,  where 
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Figure  17:  Evaluation  of  Additivity 


K 


2.5  cy/mm  . 


1 = _ J 

2 x dot  spacing  2 x .2 

Figure  17  clearly  shows  that  the  simple  sum  underestimates  the  combined 
spectrum  at  the  lower  power  levels  although  agreement  is  good  at  higher  levels. 
The  root  mean  square  deviations  from  the  combined  spectrum  were  calculated 
for  each  condition  and  are  presented  in  Table  7. 

Table  7:  RMS  Errors  in  Noise  Additivity 

(percent  of  composite  value) 

Grid  Frequency  Noise  Level  R.M.S.  Deviation 


1.16  cy/mm 

Low 

63% 

1.16  cy/mm 

High 

50% 

1 . 76  cy/mm 

Low 

55% 

1 .76  cy/mm 

High 

64% 

The  errors  are  fairly  random.  They  do  not  appear  to  be  dependent  on  noise 
level  or  grid  frequency.  The  size  of  errors  is  surprisingly  small  considering 
the  severity  of  the  conditions  imposed.  Furthermore,  the  effect  on  information 
capacity  measures  will  be  somewhat  less  since  the  major  deviations,  at  least 
in  this  case,  occur  at  the  lower  power  levels  which  in  turn  have  the  least 
effect  on  information  capacity.  The  consistency  of  the  magnitude  and 
nature  of  the  errors  (underestimation  at  low  power  levels)  also  suggest  the 
possible  development  of  an  empirical  correction.  Such  an  approach,  however, 
will  require  much  more  detailed  testing  under  more  realistic  conditions. 

Based  on  the  preceding  results,  there  is  reason  to  believe  that  for  more 
realistic  conditions  the  noise  measurement  errors  will  not  be  excessive. 

An  evaluation  like  the  one  just  discussed,  however,  should  be  performed 
for  specific  conditions  of  interest  before  such  a decision  is  made.  If 
such  an  evaluation  demonstrates  unacceptable  error  levels,  an  empirical 
correction  can  be  considered  or  an  alternate  approach  used. 
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ALTERNATE  APPROACHES 


The  information  provided  by  the  wedge  elements  in  the  ROSA  detector  offers 
some  promise  for  alternate  approaches.  The  wedges  provide  power  values  as  a 
function  of  orientation  integrated  over  a frequency  interval  determined  by 
the  detector  geometry  and  the  measurement  optics  (see  equation  (20)  in 
Section  III.  Figure  18  is  a plot  of  the  relative  wedge  values  for  the 
1.16  cy/mm  grid  with  the  high  noise  condition.  The  input  was  oriented  so 
that  the  power  spike  from  the  grid  falls  at  90°.  The  two  dimensional  noise 
pattern  produces  spikes  at  0°,  90°,  and  180°  as  well  as  on  the  diagonals, 

45°  and  135°.  The  90°  wedge,  therefore,  provides  an  integrated  spectrum 
for  the  combined  noise  and  signal.  Integrated  noise  spectra  are  provided 
by  the  wedges  at  0°,  45°,  135°  and  180°.  Unfortunately,  the  0°  and  180° 
wedges  do  not  fully  cover  the  spikes  and  the  45°  and  135°  wedges  provide 
noise  values  along  the  pattern  diagonal  and  are  only  relatively  correct  for 
the  effect  at  90°.  A 45°  tip  of  the  input  image  would  have  provided  more 
valid  data.  However  for  illustrative  purposes  signal-to-noise  values  for 
each  of  the  grid-noise  combinations  were  calculated  from 

S/N  = 90°  wedge/45°  wedge 

These  values  are  listed  in  Table  8,  along  with  the  same  calculations  from 
measurements  on  the  grid  inputs  with  noise. 


Table  8:  Wedge  Signal -to-Noise  Values 


Grid  Frequency 
(cy/mm) 

Noise  Overlay 
Level 

90°  wedge/45°  wedge 

Percent 

Loss 

1 .16 

None 

142.7 

1.16 

Low 

38.2 

73.2 

1.16 

High 

32.8 

77.0 

1.76 

None 

216.0 

1.76 

Low 

60.9 

71 .2 

1.76 

High 

59.9 

72.3 

The  percent  S/N  loss  due  to  the  overlay  is  also  shown  in  the  table.  The 
resulting  values  are  consistent  with  expectations  except  for  the  small  difference 
between  noise  levels  for  the  1.76  cy/mm  grid. 
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The  resulting  values  are  of  the  form 


/ PS,  N (kr>  7 } d kr 

/ P'N  (kr.  Y ) d kr 


and 


log2 


Jp: 5.  n (V’  d kr  f 

T e'„  <kr.  v ) d kr 


d kp  d Y 


The  wedge  signal-to-noise  values,  therefore,  cannot  be  used  directly  in  the 
information  capacity  metric  but  may  well  show  a stable  relationship  with 
capacity. 


This  approach  has  an  important  advantage  in  that  the  noise  measure  is 
internal.  The  noise  measure  is,  by  definition,  properly  representative. 

Besides  not  being  a direct  measure  of  capacity,  it  has  the  additional  disadvantage 
in  that  the  input  signal  must  be  highly  directional  (e.g.  , a line  grid). 
Furthermore,  the  present  detector  geometry  does  not  provide  a proper  range  of 
frequency  integration.  It  is  too  large,  for  example,  for  the  conditions  in 
the  test  described  above.  This  approach,  however,  remains  a potential 
technique  and  merits  further  consideration  in  Phase  II. 

A second  "indirect"  capacity  measurement  approach  using  the  wedge  data  can 
also  be  considered. 

Figure  19  compares  relative  wedge  values  for  the  1.16  cy/mm  grid  with  and 
without  the  high  noise  overlay.  It  is  intuitively  reasonable  to  consider  the 
magnitude  of  the  relationship  between  these  two  spectra  as  an  index  of  the 
"fidelity"  of  the  noi sy  response . The  correlation  between  these  spectra  (i.e., 
the  input  grid  only  and  the  output  noise  grid)  provides  such  an  index. 

Standard  linear  regression  techniques  as  outlined  in  any  statistical  textbook 
will  provide  a best  fit  equation  of  the  form 
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"s,  N ( T > 


A + B W ( T ) 


where 


“s,  N * v I 


"S  (y  > 


is  the  relative  wedge  value  at  orientation  y for 
the  output  image 

is  the  corresponding  value  for  the  input  image 


This  analysis  provides  two  additional  statistics  that  are  relevant  to  display 
performance; 

2 

The  coefficient  of  determination,  R , is  the  ratio  of  the  variance  in 

the  independent  variable  (W<-  ^ ( y ) ) that  is  accounted  for  by  the 

dependent  variable  (W<.  ( y)  ).  Conversely,  1 - R ^ represents  the 

proportion  of  the  signal  and  noise  spectrum  that  is  unrelated  to  the 
, 2 

signal.  R , therefore,  is  a form  of  signal-to-noise  ratio. 

2 

The  standard  error  of  estimate  is  derived  from  R and  is  the  absolute 
value  of  the  power  in  the  signal  and  noise  spectrum  that  is  not 
related  to  the  signal  spectrum. 

S.E.E.  = var  (WS)  N ( y ) ) (1  - R2) 

The  standard  error  of  estimate,  then,  can  be  considered  an  internal 
measure  of  the  noise  component  in  the  combined  spectrum. 


If,  for  the  grid  input,  the  90°  wedge  value  is  used  as  an  estimate  of  the 
combined  spectrum,  then  an  alternate  signal-to-noise  ratio  can  be  defined  as 


S/N 


90°  wedge 
S.E.E. 


The  Rd  and  S/N  values  were  calculated  for  this  "cross  spectral  correlation1 
approach  and  are  listed  in  Table  9. 
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Grid  Frequency 


Noise  Ovc. .'i.'.y 
Level 


1.16 
1.16 
1.76 
1.76 

? 

The  R and  S/N  values  are  me' 

agreement  in  trend.  Both  verv  . 4 ■ *: 

lower  values  for  the  1.76  cy /?..!••  •. - ir  -•  ■ .-> 

the  45°  wedge  as  a noise 
a selection  between  these  tw; 


Two  approaches  for  obtaining 
defining  and  measuring  a * 

above.  The  "off-axis  weoge ” oc; oa  • 
signal  pattern.  The  "crr.s 
or  measurement  of  the  inpu i : 
sor.e  potential  as  display  oc-r'  • 
further  consideration  ir,  Prvs- 
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SECTION  VI 

DISPLAY  PERFORMANCE  EVALUATION 


The  preceding  sections  have  dealt  with  the  basic  validity  of  a general  ap- 
proach for  the  measurement  of  display  information  capacity.  The  purpose  of 
such  a measurement  is  to  provide  a tool  for  the  evaluation  of  viewer  perform- 
ance with  imaging  displays.  A complete  description  of  performance  requires 
the  determination  of  relationships  between  levels  of  relevant  system  param- 
eters and  information  capacity.  It  is  important,  therefore,  to  identify  those 
system  parameters  that  are  likely  to  influence  information  capacity. 

The  development  of  the  basis  approach  has  identified  a number  of  display  sys- 
tem parameters  that  will  have  an  effect  on  the  optical  power  spectrum  measure- 
ments and  hence  on  the  resulting  information  capacity  value. 

SPATIAL  FREQUENCY 

Spatial  frequency  is  a principal  dimension  in  the  optical  power  spectrum  mea- 
surement. The  capability  to  determine  display  output  as  a function  of  spatial 
frequency  is  an  important  property  of  this  technique.  The  contributions,  as 
a function  of  spatial  frequency,  to  total  display  information  capacity  should 
be  a useful  discription  of  display  performance.  This  distribution  will  be 
very  sensitive  to  the  display  modulation  transfer  function. 

DYNAMIC  RANGE 

Testing  discussed  in  Section  IV  has  shown  the  information  capacity  measure  to 
be  sensitive  to  maximum  available  contrast.  The  important  parameter  here  is 
maximum  contrast  available  to  the  viewer.  Thus,  degradations  such  as  those 
caused  by  reflections  of  ambient  light  must  be  included.  The  image  capacity 
vs.  dynamic  range  distribution  will  indicate  the  sensitivity  of  a given  dis- 
play to  such  degradations. 

DISPLAY  LUMINANCE 

Average  luminance  of  the  display  will  not  have  a direct  effect  on  the  infor- 
mation capacity  measurement  because  of  zero  frequency  normalization.  However, 
the  interaction  of  luminance  with  other  parameters  such  as  dynamic  range  or 
modulation  transfer  will  be  reflected  in  the  measurement.  An  important  inter- 
action, for  example,  exists  among  all  three  parameters  for  the  CRT  displays. 
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Increased  luminance  usually  improves  the  dynamic  range  but,  since  it  also 
tends  to  increase  scanning  spot  size,  will  degrade  the  modulation  transfer 
function.  Luminance  vs.  information  capacity  should  identify  optimum  lumi- 
nance levels  for  a given  display. 

INTEGRATION  TIME 

Visual  integration  will  tend  to  reduce  the  effect  of  noise  in  dynamic  displays. 
This  factor  is  important  to  proper  measurement.  Ideally,  this  should  not  be  a 
parameter.  It  is  a characteristic  of  the  viewer  and  if  fixed  and  known  would 
be  incorporated  in  the  measurement  process.  Unfortunately,  our  knowledge  of 
visual  integrative  behavior  is  not  adequate  to  achieve  such  an  approach. 
Measurement  at  a ranqe  of  values  around  0.1  seconds  and  correlation  of  results 
with  viewer  performance  should  help  identify  a proper  value. 

CONTENT  ORIENTATION 

Directional  effects  will  occur  if  the  noise  behavior  and/or  modulation  transfer 
vary  with  direction.  This  is  likely  to  occur  with  line-scan  CRT  displays. 

It  is  generally  recommended  that  pictorial  display  resolution  be  equal  in  the 
x and  y directions.  This  recommendation  should  also  apply  to  information 
capacity.  Matrix  displays  will  exhibit  reduced  performance  on  the  diagonals. 
Information  capacity  measures  can  help  to  evaluate  the  magnitude  of  this  degra- 
dation. 

NOISE 

Noise  level  is  another  major  parameter  in  the  capacity  metric.  Noise  character- 
istics are  often  assumed  to  be  fixed  for  a given  display.  Noise  level  will,  of 
course,  vary  with  some  of  the  parameters  mentioned  above.  This  effect  should  be 
incorporated  in  the  measurement,  provided  an  appropriate  noise  sample  is  used. 

It  may,  nevertheless,  be  of  interest  to  evaluate  the  effect  of  additional,  ex- 
ternally imposed,  noise  levels  on  display  performance.  In  this  instance,  noise 
level  would  appropriately  be  used  as  a parameter  and  the  resulting  capacity  vs. 
noise  level  distribution  used  to  evaluate,  say,  the  susceptibility  to  jamming. 

INPUT  CHARACTERISTICS 

The  information  capacity  of  the  input  signal  will  influence  the  measured  display 
capacity.  Clearly,  if  the  input  level  does  not  exceed  the  capability  of  the  dis- 
play, then  display  capacity  has  not  been  measured.  Certain  input  characteristics 
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•jf-n  as  randomness  and  directionality  have  been  discussed  with  respect  to 
effect  on  capacity  measurement  and  their  utility  as  calibration  or 

‘‘•.viuremerst  Inputs. 

Measured  information  capacity  is  Greater  (at  all  frequencies)  than 
levels  of  the  information  to  be  displayed,  then  the  display  is  rated 
if  the  capacity  is  not  sufficient,  then  it  is  rated  low.  Such  a 
uegment  requires  knowledge  of  the  information  content  of  the  range  of  in- 

•j.  ■ r',  ■ the  application  of  interest. 


L CAPABILITIES 

i.e  meet  .meaningful  evaluation  must  encompass  consideration  of  both  display 
of  viewer  capabilities.  The  viewer's  performance  is  the  ultimate  criterion 
ceracT. ion  v." . fi  the  display  parameters  must  be  considered.  Visual 
time  has  already  been  discussed.  Viewing  distance  is  incorpor- 
. -i  r capacity  metric  when  spatial  frequency  is  expressed  as  cycles/ 

A thlrc;  area,  contrast  detection  thresholds,  has  been  mentioned  in 
. ..  V’. : s is  a much  more  difficult  situation.  Its  effect  can  be 

v ';-ed  by  dealing  only  with  conditions  where  display  noise  is  the  limit- 
po'-aoete’'  bu t the  problem  cannot  be  ignored  and  must  involve  a compari- 
n'" . capacity  measures  and  viewer  performance. 

v -,r  to  to1  eshc-lds,  there  is  the  problem  of  viewer  information  capa- 
riv'-e  is  iitcle  advantage  to  a high  display  information  capability  if 
:W • • cads  the  viewer.  Our  knowledge  of  visual  information  handling  capa- 
’ meager.  The  existence  of  a valid,  reliable  information  capacity 

>.r  • id  be  » useful  tool  for  extending  this  knowledge. 
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CONCLUSIONS 

a)  The  general  approach  of  optical  power  spectrum  measurement  for  toe 
determination  of  display  information  capacity  is  valid  at  ..he  oueic 

level . 

b)  The  basic  measurement  capabilities  of  the  existing  OPS  r-. 
equipment  and  calibration  procedures  are  adequate.  Rcot-mear.-squar*- 
(R.M.S.)  deviations  from  theoretical  values  are  about  10'.  R.M.S. 
measures  of  repeatability  error  are  about  2*  of  the  expected  value. 
These  error  levels  apply  over  a useful  dynamic  range  of  6 orders  or 
magnituae. 

c)  There  are  a number  of  potential  problems  in  the  practical  application 
of  the  general  approach. 

o Nonlinearity  of  amplitude  units 
o Measuring  aperture  effects 
o Phase  effects 
o CRT  film  recording 
o Noise  measurement 

d)  Promising  solutions  for  these  problems  exist,  but  most  need  further 
development  and  testing. 

e)  The  combined  influence  of  the  display  ana  viewer  characteristics  on 
information  capacity  will  include  at  least  the  following: 

o Spatial  frequency  response 
o Dynamic  range 
o Luminance 
o Integration  time 
o Input  image  characteristics 
o Noise 

o Viewing  distance 
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APPENDIX  A 
INFORMATION  THEORY 


DEFINITION  OF  IFORMATION 

A formal  definition  of  "information"  exists  within  the  framework  of  informa- 
tion theory.  This  definition  is  a precise  and  unique  measure  and  differs  in 
some  important  respects  from  common  understanding  of  information.  Informa- 
tion theory  does  not  deal  with  the  meaning  or  relative  importance  among  various 
pieces  of  information.  It  does  not  deal  with  the  formation  of  concepts  or 
decisions  by  the  operator.  The  theory  defines  information  as  a change  in 
the  state  of  knowledge.  To  illustrate,  suppose  a display  is  provided  to  in- 
dicate the  winner  in  a ten  horse  race.  If  all  horses  have  an  equal  proba- 
bility of  winning,  then  before  the  race  the  state  of  knowledge  can  be  des- 
cribed as  the  probability  of  identifying  the  winner; 

prob  (initial)  = 1/10  = 0.1 

At  the  end  of  the  race  the  display  indicates  the  winner  and,  if  we  know  that 
the  display  is  error-free,  our  state  of  knowledge  is  now  perfect,  i.e.; 

prob  (final)  = 1.0. 

The  gain  in  knowledge  or  information  is  defined  as; 

prob  (final ) 

gain  in  information  = ^ -(TnYtial)- 


For  this  display; 

gain  = = 10. 

0.1 

Suppose  a simpler  display  were  provided,  one  capable  of  indicating  only 
whether  the  number,  1 through  10,  of  the  winning  horse  was  odd  or  even. 
Then,  the  displayed  information  reduces  the  number  of  possible  winners 
to  five  and  the  resulting  state  of  knowledge, 
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prob  (final ) = 1/5  = 0.2. 


The  resulting  gain  in  information  provided  by  the  display; 


0.2 

gam  = ^ 


= 2. 


The  gain  in  information  is  equal  to  the  number  of  possible  equiprobable  display 
responses. 


prob  (final)  = M/Ny 


where,  M = number  of  possible  display  responses 
Ny  = number  of  possible  event  outcomes. 

prob  (initial ) = 1/Ny 


and 


gain 


prob  (final ) _ ^NT  _ M _ 

prob  (Initial ) 1/Ny 


1/P 


where  P = probability  of  a given  display  response. 


This  is  in  agreement  with  the  gain  values  calculated  for  the  two  displays  above. 

The  basic  unit  of  information  is  defined  as  the  “bit"  (short  for  binary  digit). 

A bit  is  the  amount  of  information  provided  by  the  answer  to  a yes-or-no  ques- 
tion when  both  answers  had  an  equal  a priori  probability  of  occurring.  This 
is  equivalent  to  the  amount  of  information  provided  by  the  odd-even  display 
above.  For  equal  probability  of  occurrence  the  information,  in  bits,  provided 
by  a given  display  response  is  defined  as; 

H(R)  = log2  (gain)  = log2  M = log2  1/p  (A.l) 

For  the  examples  above; 


73 


. 32 


b’t  3/  response 


display  H ~ loo.,  2 - 1 0 bits/response 


!')l  vi- 
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Suppose  a device  at  e -?  splay log  a series  of  responses  where  the  set 

C'T  pOS5  1 2 ! t?  :» f f; (?  j Ov  , 


V R3,. 


with  prcbaoilifits  ■■■'  :c.c or 


PiK,;,. 


Tiien  if  there 
{zero-me!Po-*y  , 
Since  ve. 
f 00,7,  a tier,  n-- 


11  correlation  between  successive  responses 

£ne  5 <f  '-air;  fry-  r:  c disp  lay  of  R.  is  log^  1/P(R. ). 

"•  " ' . oc'or--'  q is  p{R.},  the  average  in- 

: • • •'  be: 


iogo  1/p (R. } 


(A. 2) 


Shannor,  (Revs- 
probabi  ’•  i •;  .?> 


■'  S1 


" P^,)  log-  p ( R • ) 


‘ “ >"oveo  -Oos  equation  A. 2 is  maximized  when  the  response 
rsos  from  ■ovation  A.], 


(A. 3) 


na  x i iru 


i.aer  response 


/ “ 


‘ >TTi  M'i 


probabilities  are  all  equal 


NON-ZERO  MEMORY 


If  correlations  exist  between  successive  responses  a more  complex  description 
of  display  performance  must  be  used.  A Markov  model  is  appropriate  for  such 
situations.  An  mtn  order  Markov  chain  describes  a situation  wherein  the  occur- 
rence of  a response  depends  on  some  number,  m,  of  previous  responses  but  is  in- 
dependent of  those  occurring  m + 1 or  more  responses  earlier.  The  probability 
of  the  i^  response,  R^,  occurring  after  some  particular  sequence  of  m responses, 
is  written  as  p(R.j|m).  This  conditional  probability  will  have  a defined  value 
for  each  possible  previous  sequence  of  m responses,  and  similar  probabilities 
will  be  defined  for  all  the  other  R responses.  The  information  gain  will  be 
determined  by  summing  for  all  responses  over  all  previous  m sequences: 


H(R)  = 2353  p(m,  R.j ) log2  l/p(Ri|m) 
m i 

= 2><m’  Ri)  1o92  Vp(Ri|m) 
m+1 

The  last  expression  follows  because  summing  over  all  m and  i for  sequences  of 
m followed  by  R..  is  equivalent  to  summing  over  all  m + 1 sequences.  It  follows 
from  the  proof  of  expression  A. 3 that 


£p(m,  Rt)  log2  1/p (R1  Jm)  < £ p(R-)  log2  l/p(R1 ) s log2  M (A. 4) 
m+1  i 

Total  equality  holds  for  a zero-memory  sequence  with  equal  response  probabil- 
ities. The  left  hand  equality  holds  for  all  zero  memory  sequences.  Maximum 
display  information  results  when  the  display  responses  have  equal  and  inde- 
pendent probabil ities  of  occurrence  and  is  equal  to  log2  (number  of  possible 
responses) . 

To  illustrate  these  concepts,  consider  a display  that  provides  a sequence  of 
alphabetic  symbols  such  as  a teletype.  Suppose  the  display  is  operating  with 
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English  text.  The  English  language  does  not  provide  for  equal  or  independent 
probabilities  of  occurrence  of  individual  symbols.  Some  letters  occur  more 
frequently  than  others  and  the  probabilities  of  occurrence  depend  on  the  spe- 
cific sequence  of  preceding  letters.  High  order  Markov  analyses  of  the  Eng- 
lish language  are  extremely  difficult.  However,  a limiting  value  of  about 
1.6  bits  per  letter  has  been  estimated  for  the  information  (Reference  1). 

If  the  symbols  (26  letters  plus  a space)  were  not  restricted  by  the  rules  of 
English  and  occurred  with  equal  and  independent  probabilities,  then 

H(R)  = log^  M = log2  27  = 4.7  bits  per  letter. 


Thus,  a teletype  with  27  possible  responses  can  provide  information  at  a maxi- 
mum rate  of  4.7  bits  per  response.  When  the  information  is  in  the  form  of  Eng- 
lish text,  however,  the  actual  rate  is  about  1.6  bits  per  letter.  This  difference 
is  not  a function  of  the  display's  capabilities  but  is  a result  of  the  partic- 
ular information  coding  scheme  used  (i.e.,  English  text). 


{ 


CAPACITY 

For  a sequential  display  such  as  a teletype,  the  information  capacity  is  time 
related.  If  M(T)  represents  the  number  of  responses  (symbols)  displayed  in 
time  T seconds,  then  capacity  in  symbols  per  second  is 

CAPACITY  = M- 


From  expression  A. 4 maximum  information  capacity  (optimum  coding)  in  bits 
per  second  will  be 


C = 


log2  M(T) 

T 


(A. 5) 


Information  theory,  quite  properly,  expresses  the  relationship  above  as 


I 
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The  use  of  the  limit  is  to  require  that  M(T)  be  measured  over  a sufficiently 
long  time  period  to  insure  an  adequate  description  of  display  performance. 

This  is  especially  important  if  the  time  required  for  individual  symbols  is 
not  the  same  for  all  symbols  (e.g.,  Morse  Code). 

We  will  drop  the  limit  designation  in  subsequent  expressions  with  the  under- 
standing that  a proper  measurement  of  capacity  will  involve  sufficient  samp- 
ling to  represent  all  display  responses  of  interest. 

APPLICATION  TO  PICTORIAL  DISPLAYS 

Information  in  pictorial  displays  is  presented  in  a parallel  rather  than 
sequential  fashion.  The  individual  responses  are  arranged  in  a spatial  rather 
than  a time  sequence.  Furthermore,  the  spatial  arrangement  of  pictorial  dis- 
plays is  generally  two-dimensional  as  contrasted  to  the  one-dimensional  time 
sequence.  For  such  two-dimensional  spatial  cases  we  can  replace  the  time 
dimension  with  area  in  equation  A. 5. 


C = 


log2  M(D) 
D 


(A. 6) 


Capacity  is  expressed  as  bits  per  unit  area  on  the  display.  The  number  of 
different  intensity  patterns  of  area  D that  can  be  reliably  distinguished 
on  the  display  is  M(D). 


The  simplest  case  is  that  of  a matrix  of  N cells  such 
play  where  each  cell  may  be  on  or  off.  The  number  of 
terns,  M is  given  by 


as  an  LED  matrix  dis- 
different  possible  pat- 


M = 


2 


N 
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Dividing  by  the  area  of  the  matrix  we  have  as  an  expression  for  the  maximum 
capacity  of  the  display 


C = 


log2  M(D) 
D 


log2  2N  N 
D = D 


The  result  shows  that  capacity  is  simply  related  to  the  number  of  elements. 
This  is  expected  since  each  element  having  only  two  possible  responses  (on 
or  off)  provides  1 bit  of  information. 


A more  complex  situation  results  if  each  element  in  the  matrix  can  assume  any 
one  of  a number  L of  intensity  values.  Then 


M = 


L 


N 


and 


og2  N log2  L 
D 


(A. 7) 
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APPENDIX  B 

OPTICAL  POWER  SPECTRUM  MEASUREMENT 


APPROACH 

There  are  two  basic  approaches  to  the  measurement  of  the  power  spectral  den- 
sity of  an  image.  Historically,  the  most  common  approach  involves  the  mea- 
surement of  the  image  intensity  distribution,  I (x,y),  usually  with  a scanning 
microdensitometer.  The  autocorrelation  function  of  the  resulting  data  is  cal- 
culated and  Fourier  transformed  to  produce  F(k  ,k  ).  The  modulus  or  real  part 

* y 

of  the  Fourier  transform  is  then  squared  to  produce  the  power  spectrum.  This 
approach  is  time  consuming  and  expensive,  particularly  if  a two-dimensional 
spectrum  is  required.  To  meet  sampling  requirements,  an  extensive  scanning 
operation  is  required.  The  computer  requirements  (capacity  and  time)  for 
Fourier  transformation  of  a very  large  data  matrix  are  not  insignificant. 

An  alternative  approach  involves  the  use  of  coherent  optical  processing.  This 
approach  takes  advantage  of  an  important  property  of  coherent  illumination. 
Under  proper  conditions,  an  aperture  illuminated  with  coherent  light  produces 
a diffraction  pattern  that  is  related  to  the  power  spectral  density  of  the 
aperture.  With  this  approach  it  is  possible  to  sample  large  areas  rapidly 
and  to  measure  the  power  spectrum  directly  without  the  need  for  a large  com- 
puter capability.  This  approach  is  called  optical  power  spectrum  measurement 
and  is  the  one  used  in  this  study. 

FRAUNHOFER  DIFFRACTION 

The  conditions  necessary  for  the  measurement  of  optical  power  spectra  are 
those  required  for  the  production  of  a F, aunhofer  diffraction  pattern.  When 
a beam  of  light  passes  close  to  the  edge  of  an  obstacle  some  bending  of  the 
light  occurs.  All  such  departures  from  the  rectilinear  propagation  of  light 
are  called  diffraction.  The  diffraction  of  the  illuminating  plane  waves  as 
they  pass  through  a grating  is  shown  in  Figure  B-l.  The  equiphase  lines  rep- 
resent the  cylindrical  surfaces  of  the  diffracted  waves.  As  these  waves  con- 
tinue beyond  the  grating,  their  equiphase  surfaces  tend  to  blend  together  and 


look  more  and  more  like  plane  wave  fronts.  These  apparent  plane  waves  form 
moving  in  several  directions.  The  largest  plane  wave  (not  shown  here)  forms 
moving  along  the  optical  axis.  This  is  the  zero  order  wave.  The  next  largest 
waves  travel  in  directions  determined  by  the  grating  line  spacing  and  the 
light  wavelength.  The  constructive  interference  of  part  of  the  light  diffrac- 
ted down  by  the  grating  is  shown  in  the  figure.  The  path  lengths  from  consecu- 
tive grating  apertures  to  a plane  perpendicular  to  the  direction  of  the  wave 
differ  by  one  wavelength.  This  would  be  a first  order  wave.  There  is  also  a 
similar  wave  diffracted  up.  Second  order  waves  are  formed  by  the  construc- 
tive interference  of  light  with  aperture  path  lengths  differing  by  two  wave- 
lengths. 

As  shown  in  Figure  B-l,  the  angle  of  diffraction,  a,  is  related  to  the  wave- 
length, X,  and  the  line  frequency,  k,  by  the  formula 

Sin  a = Xk  (B.  1 ) 


for 


- n/2  < a < Tt /2 

A plane  wave  diffracted  by  a line  grating  forms  a series  of  plane  waves  mov- 
ing off  at  angles  predicted  by  the  above  formula.  As  shown  in  Figure  B-2, 
a plane  wave  passing  through  a lens  is  focused  to  a spot  in  a plane  located 
one  focal  length  behind  the  lens.  The  diffraction  pattern  formed  in  this  plane 
is  called  the  Fraunhofer  diffraction  pattern.  The  position  (r)  of  the  spot 
in  the  focal  plane  (Fraunhofer  Diffraction  Plane)  is  taken  into  account  by  the 
equation 

tan  a = -p- 

where  F is  the  focal  length  of  the  lens. 

Since,  from  equation  (B.l) 


Sin  a = Xi< 


80 


and 


sin  a s tan  a for  | a | < 10°  (B.2) 

r = Fx  k 

Therefore  a spatial  frequency  k in  the  input  image  will  cause  light  to  be  con- 
centrated in  the  diffraction  plane  at  a radius  equal  to  Fxk.  The  output  posi- 
tion is  only  a function  of  the  diffracted  angle.  It  is  independent  cf  x,  y 
position  in  the  input  plane  but  is  sensitive  to  orientation. 

THE  DIFFRACTION  INTEGRAL 

The  light  amplitude  at  any  point  in  the  Fraunhofer  diffraction  pattern  is  pro- 
portional to  the  Fourier  coefficient  for  the  corresponding  spatial  frequency 
component. 

Specifical ly, 

A(x,y)  exp  £-2tti7x(six  + my)  J dx  dy  (B.3) 

where  1 and  m are  angular  coordinates  defined  by  the  diffraction  angles 

(p.  = sin  «x;  m = sin  a ) 

A(t,m)  is  the  diffraction  pattern  amplitude  distribution 
x is  the  wavelength  of  the  coherent  illumination 
F is  the  focal  length  of  the  transform  lens. 

The  limits  of  integration  are  defined  by  the  input  image  area.  The  derivation 
of  equation  (B.3)  is  somewhat  lengthy  and  will  not  be  treated  here.  The  topic 
is  covered  in  most  standard  optical  texts.  References  2 and  11  are  recom- 
mended. 

A(x,y)  is  the  amplitude  distribution  in  the  input  plane.  It  can  be  expressed 
as 

A(x,y)  = LAT(x,.y)  (B.4) 


A(?,m)  aF  JJ 


82 


f 


‘ 'J'J  » m iu  i 


where  L = the  amplitude  of  the  coherent  illumination  (assumed  to  be  constant 
over  the  image  area) 

A-p(x,y)  = amplitude  transmittance  distribution  at  the  input  image  plane. 


So, 


/1(c,m)  = 


Lp  If 


AT(x,y)  exp  | -2tt/a  (f.x  + my)  dx  dy  (B.5) 


F rom  B . 1 


i = sin  «x  = Akx,  and 
m = s i n » = A ky 


and 


A(rx,ry ) = fp  j y"  AT(x>.y)  exp  -2ni(kxx  + kyy)J  dx  dY 


where  r and  r are  linear  dimensions  in  the  diffraction  plane, 

* y 


r = {/ F,  and 
x 

ry  = m/F. 


The  integral  is  the  Fourier  transform  of  Ay(x,y)  so. 


A^rx,ry)  ' aF  FT  (kx,ky)‘ 


(B.  6) 


(B.  7) 
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The  amplitude  distribution  of  the  diffraction  pattern,  however,  cannot  be 
measured  directly.  Photodetectors  measure  intensity,  I,  rather  than  ampli- 
tude. Since 
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